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AFOSI 


The  Martian  and  Venusitn  m^^el  »t|Tiosphcres  utilized 


followed  those  given  in  a  p9elimiAtL9;|  engineering  estimate^ 

As  the  charactefistics  of  these  atmotplte#es  are  not  known  precisely,’  probable 
bounds  were  established.  These  bo»n<*s  for  the  Martian^atmospherg^re 


t  abi-ein  >^4  w.t-rr  'Tabi^ 


V>)cU. 


euntaiwn^hc  planetafj|  d^^fr^c  (l4i|a  the  two  planet^#  As  the  'Venusian 
^Xtmos^erc  dlTta  was  ^ive#  4#  teitu*  ©f  temperature  and  pressure, 


the  density  was  calcwlaneii 


and  is  plotted  as  a  function  of 
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TABLE  1  - 1 .  MARS  ATMOSPHERE* 
INDEPENDENTLY  DETERMINED  PARAMETERS 


Symbol  a 

Parameter 

Unit 

Low 

Nominal 

High 

So 

Acceleration  of 

ft/aec^ 

12 

12 

13 

Surface  Gravity 

cm/aec^ 

360 

375 

390 

•^m 

Surface  Maaa 

slug /ft  2 

9.  5 

16 

22 

Preaaure 

gm/cm 

150 

250 

350 

T 

Stratoaphere 

OR 

230 

320 

410 

Temperature 

OK 

130 

180 

230 

To 

Surface 

OR 

*  380 

410 

470 

Temjlerature 

OK 

210 

230 

260 

Cp 

Specific  Heat 

BTU/ 

7.  4 

7.7 

8.  1 

at  Constant 

slug  OR 

Preaaure 

cal/gm  Ok 

.23 

.24 

.25 

♦Tentative 

• 

'  • 

* 

TABLE  1-2.  MARS  ATMOSPHERE 

parameters  dependent*  on  selected  composition  of  ATMOSPHERE 

Symbol  Parameter  Unit  Low  ^  ■  (g)  — — ^  High 

High  ^  —  (Mm)  — ^  Low 

Cg  Concentration  of  % 

Conatituenta  by 
Volume 


COg 

0.7 

1.9 

7.  2 

A 

0.6 

6.  0 

©2**  . 

0 

-- 

Nz 

98.  7 

86,8 

Concentration  of 
Conatituenta  by 

Maaa 

% 

CO2 

1.  1 

3 

11 

A 

0.8 

8.7 
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TABLE  1-2.  MARS  ATMOSPHERE  (Conf  d) 

PARAMETERS  DEPENDENT*  ON  SELECTED  COMPOSITION  OF  ATMOSPHERE 


Symbol 

Parameter  Unit 

Low  ^ 
High  V 

0 

-  (g)  • 

-  (Mjal 

- ^  High 

M  Low 

(2.7) 

N2 

98.  1 

80.  3 

M 

Average  Molecular 

Weight  by  Volume 

28.  2 

28.  4 

29.9 

y 

Ratio  of 

1.  42 

1.  41 

1.  41 

Specific  Heats 

High  ^ 

(g»  ^m) 

- )  Low 

P 

Surface  Pressure  Ib/ft^ 

284 

196 

113 

mb. 

136 

94 

54 

’i’Data  are  consistent  within  columns 

♦♦Conjectural  Oxygen  not  considered  in  Average  Molecular  Weight 


TABLE  1-3.  MARS  ATMOSPHERE 
DERIVED*  PARAMETERS 


Symbol 

Parameter 

Unit 

Lowf_g, 

Po,  M,  “ 

— >  High 

Highf-T 

— ^  Low 

^o 

Surface  Density 

slug  X  10"^ 

1.4 

2.7 

4.  3 

ft^ 

gm  X  10*'^ 

0.  75 

1.4 

2.  2 

cm^ 

H 

Scale  Height 

ft  X  10^ 

59 

41 

29 

Km 

18 

12.  5 

9 

Reciprocal  of  H 

1  X  10-5 

1.7 

2.4 

3.  4 

■"TT 

1  X  10-2 

5.6 

8.0 

11.  0 

Km 

T 

Dry  Adiabatic 

OR  X  10-5 

-1.96 

-2.04 

-2.  13 

Lapse  Rate 

ft 

OK 

-3.  43 

-3.75 

-4.  10 

Km 
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TABLE  1-3,  MARS  ATMOSPHERE  (Cont'd) 
DERIVED*  PARAMETERS 


Symbol 

Parameter 

Unit 

Low  ^ —  g, 

Po,  M| 

4  High 

So 

Speed  of  Sound 

ft/ sec 

High  i — T_ 

1064 

1013 

Low 

985 

at  Surface 

m/  sec 

325 

309 

300 

8 

Speed  of  Sound 

ft/  sec 

1013 

903 

765 

u 

in  Stratosphere 

m/  sec 

309 

275 

233 

T 

Lapse  Rate 

1  X  10’^ 

-4.  17 

-4.  98 

-5.  61 

T 

^o 

Surface  Temp. 

ft 

1  X  10-2 

-1.  32 

-1. 63 

-1.  95 

Km 

*  Data  are  consistent  within  columns 
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Figure  1.  1  Temperature  of  Mare  Atmosphere, 


Figure  1.2  Density  of  Mars  Atmosphere. 
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TABLE  1-4.  COMPOSITION  OF  VENUSIAN  ATMOSPHERE 


I  ("best") 

II 

III 

90%  N2 

96%  N2 

74%  N2 

9%  CO2 

3%  CO2 

25%  CO2 

1%  A 

1%A 

1%  A 

m  s  29. 6 

m  s  28.  ^ 

m  B  32. 1 

5  =  1.39 

6  =  1.40 

6  =1.38 

A  ■ 

JL 

■A '  ’ 

T  =  -8,  5®K/km 

r=  -8.4®K/km 

r  =  -8.9®K/km 

m  s  mean  molecular  mats 
6  r  mean  ratio  of  specific  heats** 
r  s  dry  adiabatic  lapse  rate* 


^These  quantities  are  given  for  STP.  At  the  temperatures  and  pressures 
which  appear  to  be  present  near  the  surface  of  Venus  they  are  badly  in 
error. 
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7 
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300.0  9.65  x  10’^®  1.457  x  10-25  150  29.6  1.40  796 


TABLE  1-8.  PLANETARY  DYNAMIC  DATA  --  VENUS  AND  MARS 


Venus 

Mars 

Radius 

0.  97  +  0.01  Rg 

6200  +  100  Km 

0.  52  R^ 

3400  +  50  Km 

Mass 

0.82  +  0.01  Mg 

4.  90  X  1024  Kg 

0.  107  Mg 

6.40  X  10^^  Kg 

Density 

5.  0  +  0.  3  gm/cc 

4,  0  +  1  gm/cc 

Oblateness 

1% 

1/192 

Surface  Gravity 

0.  88  +  0.  03  g  2 
860  +  30  cm/  sec 

390  +  10  cm/ sec 
0.  39'g 

Orbit  Details 

Mean  distance  to  Sun 

108.  1  X  10^  Km 

227.  8  X  10^  Km 

Maximum  distance  to  Sun 

109.4  X  10^  Km 

248  X  10^  Km 

Minimum  distance  to  Sun 

106.  5x10^  Km 

206  X  10^  Km 

Eccentricity  of  Orbit 

0.  007 

0.093 

Inclination  of  orbit  to  plane 
of  ecliptic 

3°  24' 

1°  51' 

Minimum  distance  to  Earth 

42  X  10^  Km 

1°  51' 

Maximum  distance  to  Earth 

257  X  10^  Km 

Sidereal  period 

225  days 

687  days 

Synodic  period 

584  days 

780  days 

Mean  orbit  velocity 

35.  0  Km/  sec 

24.  2  Km/  sec 

Semi -major  axis 

108.  21  X  10^  Km 

227.9  X  10^  Km 

period  of  rotation 
escape  velocity 
inclination  of  Equator 
to  orbital  plane 

Satellites 

undetermined 

10.  4  Km/  sec 

undetermined 

none  detected 

24  hrs  37  min 

5.  0  Km/  sec 

25°  12' 

2 
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DENSITY,  SLUGS /FT 


Figure  1,  3  Venusian  Model  Atmospheres. 
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II.  PROGRAM  OBJECTIVES 


The  objective  of  this  comprehensive  program  of  research  was  to  solve 
the  problem  of  ballistic  entry  into  the  atmosphere  of  Mars  and  to  establish 
the  techniques  for  its  solution  for  Venusian  entry.  This  requires  the  theo¬ 
retical  and  experimental  determination  of  the  heat  transfer,  the  aerodyna¬ 
mic  forces, and  the  behavior  of  typical  thermal  protection  systems  during 
the  ballistic  entry  of  a  space  probe  into  the  atmosphere  of  Mars. 

The  investigators  in  the  Electron  Physics  Section  performed 
studies  in  order  to; 

1.  Prepare  shockwave  data  and  Mollier  tables  for  Mars  and  Venus. 

2.  Experimentally  determine  dissociation  rates  typical  for  the  atmo¬ 
sphere  of  Mars. 

3.  Construct  and  investigate  the  feasibility  of  using  an  arc -heated 
helium  driven  shock  tube-tunnel  for  the  determination  of  heat  transfer  and 
aerodynamic  data  by  evaluating  and  calibrating  a  prototype  facility,  and 
developing  an  adequately  sized  test  facility  utilizing  the  prototype  results. 

4.  Determine  analytically  the  shock  facility' s  expected  performance 
with  N2  and  CO2. 

5.  Develop  new  transducers  and  instrumentation  concepts  for  acquir¬ 
ing  model  data  in  the  high  enthalpy  shock  facility  (particularly  for  obtaining 
heat  transfer). 

6.  Determine  experimentally  the  aerodynamic  heat  transfer  to  the 
stagnation  point  of  a  blunt  body  during  ballistic  entry  into  the  atmosphere 


12 


of  Mars  and  investigate  the  heat  transfer  for  similar  entry  into  Venus. 

7.  Establish  the  beha'/ior  of  various  ablation  shield  materials  in  en¬ 
vironments  typical  of  hypervelocity  entry  into  atmosphere  of  Mars. 

8.  Determine  Venusian  ionosphere  attenuation  effects. 


III.  STUDY  RESULTS 


III-l  EQUILIBRIUM  THERMOCHEMISTRY* 

For  both  the  theoretical  and  experimental  studies  of  planetary  entry 
it  was  necessary  first  to  establish  the  conditions  behind  the  shock  wave. 

This  required  the  determination  of  the  equilibrium  gas  composition.  Two 
approaches  were  taken  to  the  problem:  first,  an  older  method  made  use  of 
equilibrium  constants  to  express  abundance  of  certain  chosen  species  in 
terms  of  the  abundance  of  other  aribtrarily  chosen  (or  more  abundant)  specie 
This  method,  referred  to  as  the  Brinkly  method  (Ref.  3),  depends  entirely 

on  the  solution  of  chemical  equilibrium  equations.  As  a  consequence,  the 

• 

program  for  the  Brinkly  method  must  be  readjusted  for  each  different  set 
of  elemental  substances  in  a  mixture,  i.e.  ,  a  (Brinkly)  program  written 
for  an  N,  O  system  is  not  applicable  to  an  H,  O  system.  However,  as  the 
program  for  such  computations  was  nearly  complete  and  checked  out  for  a 
number  of  species  and  mixtures,  it  was  utilized  to  obtain  early  results. 

This  work  supported  by  JPL  contract  was  completed  (Refs.  4  and  5)  and 
was  used  throughout  other  portions  of  this  and  other  studies.  A  second 
approach  to  this  problem  was  undertaken  as  a  part  of  this  study.  This 
approach,  referred  to  as  the  Rand  method,  makes  no  distinction  among  con¬ 
stituent  species.  It  is  based  on  the  fact  that  any  n  component  chemical  sys¬ 
tem  is  in  thermodynamic  equilibrium  when  the  free  energy  of  that  system 
(i.e.  ,  the  Gibbs  function)  is  at  a  minimum.  With  this  procedure  where  the 
free  energy  equation  is  linearized,  the  input  required  consists  of  values  of 

"iPrincipal  Investigators:  M.  J.  Linevsky  and  H,  Y.  Wachman 
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the  standard  molar  free  energy  of  each  constituent  and  the  quantities  of  each 
elemental  substance  in  the  mixture.  With  this  approach,  four -decimal  place 
accuracy  in  constituent  mole  fractions  can  be  obtained.  Mole  fractions  of 
trace  constituents  can  be  calculated  by  using  a  simple  subroutine  which  im¬ 
proves  the  computed  mole  fractions  by  substituting  them  in  a  set  of  linearized 
chemical  equilibrium  equations.  In  this  case,  the  chemical  equilibrium 
equations  are  only  used  to  improve  the  mole  fractions  and  obtain  values  for 
minor  constituents.  If  the  elemental  substances  of  the  mixture  change,  the 
procedure  requires  only  a  change  in  the  data  input  and  in  the  afore  mentioned 
trace  element  subroutine.  The  latter  change  is  not  required  if  only  major 
constituents  and  thermodynamic  properties  are  desired.  Simple  algebraic 
forms  suffice  for  expressing  the  free  energy  of  constituent  or  mixture,  and 
a  variety  of  computational  methods  are  available. 

In  addition  to  the  linear  programming  method,  the  free  energy  equa¬ 
tion  can  be  minimized  using  a  desk  computer  by  a  method  of  steepest  de¬ 
scent.  This  has  already  been  worked  out  in  a  way  consistent  with  the  linear 
programming  approach. 

In  general  the  Rand  method  can  be  readily  adapted  to  constant  volume 
and  constant  pressure  problems;  the  former  which  is  usually  quite  complex 
leads  to  an  even  simpler  formulation  in  this  case  than  the  latter. 

Extention  of  the  method  to  multiphase  systems  and  to  ionic  systems 
was  straightforward. 

Trace  elements  can  be  eliminated  or  included  in  computations  at  will 
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and  this  does  not  hamper  the  accuracy  of  the  final  result.  Even  if  a  trace 
element  has  been  omitted  in  the  original  formulation,  it  can  be  included  by 
simple  calculation. 

The  method  was  also  modified  so  that  the  mole  fraction  of  any  consti¬ 
tuent  can  be  fixed  in  advance,  and  the  equilibrium  composition  and  thermo¬ 
dynamic  properties  of  a  multicomponent  system  containing  a  fixed  mole 
fraction  of  any  constituent  are  calculable.  For  example,  the  method  can  be 
used  to  compute  the  properties  of  a  chemical  system  in  which  a  certain 
ionization  level  is  to  be  achieved  by  seeding  with  a  low  ionization  potential 
material.  The  desired  ionization  level  is  introduced  as  an  input,  and  the 
program  is  used  to  calculate  the  equilibrium  composition  and  the  quantity 
of  material  required  to  achieve  equilibrium  In  general,  convergence  is 
quite  rapid  regardless  of  starting  point  and  thus  requires  little  machine 
time.  The  computer  program  has  been  completed  for  a  multicomponent 
chemical  system  including  ionized  species  and  has  been  augmented  to  take 
advantage  of  the  maximum  capacity  of  the  IBM  7090  computer.  The  sub¬ 
routing  for  trace  constituents  was  written  based  on  a  list  of  components 
representative  of  planetary  atmospheres.  Subsequent  to  a  check  out,  this 
program  will  be  available  for  fast  inexpensive  calculations  of  chemical 
equilibrium  compositions  and  thermodynamic  properties  for  any  perturba¬ 
tion  of  the  still  very  uncertain  planetary  atmospheres  or  test  gases  used  to 
represent  the  planetary  atmosphere  in  the  ground  simulation  test  facilities. 
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ni-2  DISSOCIATION  RATE  STUDIES* 


Studies  of  non-equilibrium  processes  in  various  gases,  especially 
nitrogen,  have  been  carried  out  over  the  past  two  yearsnn  a  combustion 
driven  stainless  steel  2"  x  2"  shock  tube  in  which  good  purity  control  is 
maintained  (impurities  ^20  ppm).  Figure  2.  1  shows  the  tube  and  associated 
equipment.  A  Mach-Zehnder  interferometer  capable  of  0.  1  micro  second 
time  resolution  was  used  to  follow  the  changing  density  in  the  dissociation 
zone.  Knowledge  of  the  shock  strength  together  with  the  hydrodynamic 
constraint  equations  plus  an  assumed  kinetic  mechanism  led  to  inferred 
values  for  the  dissociation  rate  constant.  An  f-3  time-resolved  spectro¬ 
graph  was  used  to  simultaneously  examine  and  record  the  white  and  mono¬ 
chromatic  fringes  of  the  emitted  light  from  the  reaction  zone.  Succt‘  iful 
control  of  the  combustion  technique  allowed  observation  of  high  degree  <  of 
dissociation  (up  to  70%)  in  relatively  rich  mixtures  at  relatively  high  ini  al 
pressures. 

The  past  work  accomplished  with  this  facility  has  been  documented  in 
reports  to  the  sponsoring  agencies;  Ballistic  Systems  Division  and  Office 
of  Scientific  Research  of  the  U.  S.  A.F.  Measurements  of  the  types  which 
have  been  conducted  are  required  to  assess  the  importance  of  non-equili¬ 
brium  processes  in  planetary  atmospheres. 

In  the  nitrogen  study  (Ref.  6)  two  points  were  stressed.  The  first 
point  emphasized  was  an  acceptable  degree  of  purity.  By  actual  mass 
spectrometric  analysis  of  gas  samples  taken  from  the  tube  it  was  found  that 
*  Principal  Investigator:  B.  B.  Cary 
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impurity  levels  were  leas  than  20  ppm.  The  second  emphasis  was  placed 
on  demonstrating  that  steady,  uniform  flow  fields  could  be  generated  by 
combustion  driven  shocks,  providing  great  care  was  taken  to  match  the 
peak  burning  pressure  of  the  combustible  load  to  the  burst  pressure  of  the 
diaphragm.  Figure  2.  2  shows  a  white  light  interferogram  taken  of  a  Mach 
10.  3  shock  propagating  into  a  13%  N2  -  87%  Ne  mixture.  The  shock  arrives 
at  the  left  at  the  upstream  interferometer  beam  causing  an  upward  jump  of 
the  fringes  followed  by  a  further  rise  in  the  dissociation  zone.  The  fringes 
then  show  only  a  small  decay  until  the  time  when  the  shock  arrives  at  the 
downstream  beam.  Here  the  fringes  revert  to  their  original  position  and 
show  a  good  bit  of  mirror  symmetry  in  the  dissociation  zone.  This  sym¬ 
metry  attests  to  the  control  of  the  combustion  technique.  It  should  be 
pointed  out  that  a  monochromatic  interferogram  can  be  simultaneously 
recorded  with  the  white  light  picture.  Figure  2.  3  shows  a  time-resolved 
spectrogram  of  the  luminosity  emitted  behind  a  M  12  shock  in  a  mixture  of 
1%  H2  -  99%  Ar.  The  stark  broadening  of  the  H/5  line  is  the  most  promi¬ 
nent  line  in  the  center.  It  is  considerably  broadened  owing  to  the  high  ion 
densities  attained  (10^^  per  cc).  All  other  visible  lines  have  been  identified 
as  argon  lines. 

Ten  theoretical  density  profiles  have  been  generated  by  a  program 
written  for  the  IBM  1620.  This  program  allows  the  computation  of  the 
density  profiles  behind  normal  shocks  resulting  from  the  bimolecular  dis¬ 
sociation  of  a  diatomic  gas  in  a  rare  gas  mixture.  The  experimentally 
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determined  profiles  were  fitted  to  the  theoretical  profiles  in  order  to  find 
the  collision  efficiencies  of  and  N  in  promoting  N2  dissociation.  Addi¬ 
tional  iteration  was  undertaken  to  obtain  final  values.  Experimentally,  a 
perturbation  on  the  N2  and  N  system  will  be  made  later  by  operating  with 
nitrogen  and  small  amounts  of  C02>.  These  results  will  be  compared  to  the 
nitrogen  data  to  determine  the  influence  of  the  presence  of  the  small  quan¬ 
tities  of  CO2  es  in  the  atmosphere  of  Mars. 

The  final  expression  for  the  rates  of  nitrogen  dissociation  as  deter¬ 
mined  in  the  recent  study  are; 

N2-N?:  k^  =  (6.  8  X  10^^)  X  T‘^*^  X  exp  +  25% _ SS. _ 

(N2-N):  k^  =  (5.  5  x  10^9)  X  T"^-®  X  exp  +  40%' 


cc 


mol-sec. 


(Np-Ne):  kj  =  (1.  5  x  10^^)  x  T"^*  ®  x  exp  +  40%  - — 

*  “3  T  ■  mol-i 


T 

-113,260 


mol-sec. 


(N2-A):  k^^  =  (2.  1  X  10^^)  X  T"^*  ®  X  exp  (  ^ 


)  +  40% 


cc 


mol-sec. 

At  6,  400°K  the  nitrogen  atom  is  about  as  efficient  as  the  nitrogen 
molecule  in  promoting  dissociation.  The  vibrational  and  dissociative  zones 
were  uncoupled  over  the  temperature  and  density  ranges  covered.  Photo¬ 
graphing  the  luminosity  behind  M  14  shocks  in  pure  nitrogen  with  an  f-3 
spectrograph  failed  to  reveal  the  1  ~  band  system  of  N2^. 

Table  3.  1  gives  a  comparison  of  the  inferred  recombination  rates  of 
this  study  with  previous  results  by  other  investigators. 
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Table  3-1 


INFRARED  RECOMBINATION  RATES  ^ _ 

2 

mole  -  sec 


Source 

Third  Body 

300°K 

6,  000°K 

This  paper 

1.  6 

4.  2 

10^“^ 

N 

2.0 

lo'^ 

4.7 

lo'^ 

Ne 

5.4 

10^5 

1.  3 

10^4 

A 

8.  0 

10^5 

1.9 

10^4 

7 

Byron 

^2 

8.  6 

10^^ 

4.4 

10^4 

N 

1.9 

10*6 

2.  0 

10^5 

A 

1.  5 

10^® 

1.  6 

10^4 

Eckerman 

Xe 

5.  1 

10^^ 

2.7 

10^5 

9 

Camm 

^2 

1.  1 

10^5 

N 

.... 

8.7 

io'5 

Allen^° 

N 

6.  5 

10^5 

Herron 

^2 

5.7  10 

(2970K) 

A 

3,  3  10 

(297°K) 

,  12 

Wentinck 

Nz 

1.  2 

ioi«> 

Harteck^  ^ 

N2 

6.  1 

10^5 

A 

6.  1 

10^^ 

»  *  •  « 

20 


21 


Figure  2.  1  Shock  Tube  end  Associated  Instrumentation 
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Figure  2.  2  White  Light  Interfere  gram  of  a  M  10.  3  Shock  Wave. 
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Figure  2.  3  Time  Resolved  Spectrogram. 
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ni-3  EXPERIMENTAL  HEAT  TRANSFER  STUDIES* 


3.  1  Introduction 

One  of  the  many  problems  connected  with  space  flight  is  the  intense 
heating  experienced  by  the  vehicle  during  its  entry  into  the  planetary  atmo¬ 
sphere.  During  the  flight  at  relatively  low  altitudes  a  shockwave  forms 
ahead  of  the  body.  The  free  stream  gas  passes  through  this  shock  wave 
and  undergoes  a  large  increase  in  internal  energy  which  is  sufficient  to 
cause  dissociation  of  molecules  and  even  partial  ionization  of  the  gas.  This 
hot,  chemically  reacting  gas  flows  around  the  vehicle  and  heat  is  transferred 
to  its  surface.  Two  modes  of  heat  transfer  can  be  distinguished.  One  is 
the  aerodynamic  heat  transfer  process  which  includes  conduction,  convec¬ 
tion  and  diffusion  of  molecules,  atoms,  ions  and  electrons  across  the  bound¬ 
ary  layer;  the  other  is  the  radiation  from  the  hot  gas  cap  (shock  layer) 
surrounding  the  forward  section  of  the  vehicle. 

The  purpose  of  this  experimental  study  was  to  investigate  the  aero¬ 
dynamic  heat  transfer  to  the  stagnation  point  of  a  blunt  body  in  a  hypersonic 
flow  of  simulated  atmospheres  of  Venus.  Besides  supporting  the  available 
analytical  studies,  the  results  of  this  work  provide  heat  transfer  data  for 
a  range  of  flight  conditions  beyond  the  scope  of  the  present  analyses. 

It  can  be  shown  that  at  the  stagnation  point  of  a  blunt  body  in  hyper¬ 
sonic  flight  the  aerodynamic  heating  is  accurately  simulated  if  the  condi¬ 
tions  at  the  edge  of  the  boundary  layer  including  the  pressure  distribution 
and  wall  conditions  are  correctly  reproduced. 

’’Principal  Investigators:  J.  S.  Gruszczynski,  C.  J.  Harris  and  W.  R. 
Warren 
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The  experiments  were  carried  out  in  the  newly  developed  arc -heated 
2"  shock  tube  which  is  capable  of  producing  in  the  laboratory  the  environ¬ 
mental  conditions  encountered  by  the  entry  vehicle.  The  correct  conditions 
of  enthalpy  and  stagnation  pressure  were  simulated  in  two  assumed  compo¬ 
sition  models  of  Venusian  atmosphere  for  flight  velocities  up  to  about  40,  000 
ft/ sec. 

Preliminary  analytical  work  involved  the  calculations  of  shock  tube 
performs. ._e  for  both  test  gas  mixtures  corresponding  to  the  specified  at¬ 
mosphere  of  Venus.  Thermodynamic  properties  behind  incident  and  normal 
shocks,  and  at  the  stagnation  point  of  the  hemispherical  body  were  computed. 
An  estimation  of  the  radiation  heat  flux  was  also  made  to  assess  the  magni¬ 
tude  of  radiant  heating  sensed  by  the  stagnation  point  heat  gage.  To  con¬ 
firm  theoretical  calculations  an  attempt  has  been  made  to  measure  the 
radiant  heating  using  a  new  type  of  a  total  radiation  gage.  The  primary 
objective  of  these  radiation  measurements  was  to  demonstrate  that  the  radi¬ 
ative  contribution  to  the  measured  heating  results  was  small,  that  is,  they 
are  used  to  verify  that  the  measured  heat  transfer  rates  are  essentially 
caused  by  the  convective  or  boundary  layer  heat  transfer  mechanism. 

3.  2  Test  Facility 

The  shock  tube  has  been  used  for  many  years  in  experiments  for  the 
generation  of  aerodynamic  data  in  flows  of  high  energy  and  high  stagnation 
pressure  (Refs.  14,  15  and  16).  The  demand  for  an  experimental  facility 
in  the  area  of  simulation  of  orbital  or  even  super-orbital  enthalpies  and 
pressures  forced  further  development  of  its  performance.  It  was  obvious 
that  this  was  beyond  the  capability  of  a  conventional  combustion  driven  shock 
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tube.  The  arc-heated  shock  tube  used  for  the  present  experiments  was  de¬ 
veloped  specifically  for  that  purpose  (Ref.  17).  It  was  originally  designed 
as  a  prototype  of  a  large  facility  but  has  been  used  successfully  in  experi¬ 
mental  investigations.  Numerous  heat  transfer  tests  were  run  with  air, 

CO2  and  as  the  test  gas  and  some  of  the  data  were  reported  in  Reference 
16. 

Figure  3.  1  shows  schematically  the  design  details  of  the  whole  shock 
tube  assembly.  The  driven  section  is  2  in.  in  diameter  and  it  is  135.  5  in. 
long.  The  driver  is  17  in.  inside  length  with  2  1/4  in.  internal  diameter. 

A  Lexan  cylindrical  sleeve  with  1/4  in.  wall  thickness  fits  inside  the  driver 
tube.  It  is  used  to  provide  electrical  insulation  along  the  arc  path;  that  is, 
it  insures  that  the  arc  discharge  is  contained  within  the  driver  gas  and  that 
the  energy  is  uniformly  distributed  axially  along  the  driver.  A  thin  wire  is 
used  to  initiate  the  arc  discharge.  This  is  located  on  axis  so  that  the  wall 
losses  will  be  uniform  in  the  axial  direction. 

The  high  voltage  electrode  which  is  insulated  from  the  driver  end 
flange  by  a  Lexan  bushing  is  in  the  form  of  a  hollow  copper  cylinder  with  a 
3/8  in.  diameter.  A  1/8  in.  diameter  brass  rod  insulated  by  teflon  tubing 
and  held  in  position  by  a  tapered  steel  washer  and  a  nut  is  inserted  through 
the  center  of  the  main  electrode.  A  5  mil  tungsten  wire  is  attached  to  the 
spider  electrode  at  the  diaphragm  end  of  the  driver  and  by  about  1  in.  cotton 
thread  to  the  copper  main  electrode.  The  position  of  the  brass  rod  relative 
to  the  main  electrode  determines  the  breakdown  voltage.  This  adjustment 
is  made  with  the  main  capacitors  bank  disconnected,  using  a  small  capaci¬ 
tor  to  prevent  breaking  of  the  wire. 
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The  electrical  circuit  for  the  driver  heating  system  is  shown  in  Figure 
3.  2. a.  The  energy  for  the  arc  discharge  is  stored  in  forty-seven  capacitors, 
each  rated  at  5  microfarads  and  20,  000  volts. 

Figure  3.  2.b  shows  a  typical  current  wave  form.  It  indicates  that  the 
circuit  is  critically  damped  and  that  the  energy  is  completely  dissipated 
in  about  50  microseconds,  and  before  the  diaphragm  opens  which  occurs 
about  150  microseconds  later.  Full  description  of  the  considerations  enter¬ 
ing  into  the  design  of  this  test  facility  is  presented  by  Warren,  et  al  (Ref.  17). 

Helium  is  used  as  the  gas  heated  in  the  driver.  The  diaphragm  sep¬ 
arating  the  driver  and  the  driven  section  of  the  shock  tube  is  made  of  stain¬ 
less  steel.  It  is  scribed  to  a  controlled  depth  (generally  20%  -  40%  of  the 
thickness)  along  two  diametrical  cross  lines.  The  required  conditions  of  a 
given  test  run  dictate  the  proper  selection  of  the  diaphragm.  The  diaphragm 
is  found  to  open  cleanly  with  a  negligible  loss  of  material. 

As  indicated  in  Ref.  17  ,  an  extensive  research  program  was  neces¬ 
sary  before  a  satisfactory  operation  of  this  facility  was  achieved.  A  photo¬ 
graph  of  the  complete  test  facility  including  the  power  supply,  vacuum  pumps, 
etc,  ,  is  shown  in  Figure  3.  3.  The  performance  of  the  tube  in  air  is  shown 
in  Figure  3.  4  where  the  attainable  shock  velocity  is  plotted  as  a  function  of 
initial  pressure  and  energy  input.  The  experimental  data  indicate  that  shock 
velocity  up  to  30,  000  ft/ sec  can  be  obtained  with  1  mm  Hg  driven  tube  initial 
pressure  and  full  energy  in  the  capacitors.  With  a  lower  initial  pressure, 
shock  velocity  up  to  36,  000  ft/ sec  can  be  reached.  Superimposed  is  a  curve 
depicting  the  performance  of  a  combustion  driver  facility.  The  limitation 
of  the  latter  is  self-evident. 
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3.  3  Experimental  Technique 


Large  amounts  of  carbon-like  residue  are  normally  found  to  cover 
the  driver  tube  walls  and  models  after  each  test.  The  origin  of  this  is  the 
Lexan  sleeve  used  to  insulate  the  driver.  It  appears  that  during  the  electri¬ 
cal  discharge  the  inner  wall  surface  of  the  sleeve  decomposes  slightly  and 
particles  of  carbon  so  formed  are  carried  along  the  driven  tube  by  the  ex¬ 
panding  driver  gas. 

To  insure  a  low  degree  of  contamination  of  the  test  gas  several  steps 
were  taken  after  each  test  run  and  before  a  new  run  was  made.  The  driven 
tube  was  thoroughly  wiped  with  alcohol  until  no  trace  of  carbon  could  be  no¬ 
ticed  on  the  inner  walls.  The  same  applies  to  the  test  model,  which  was 
also  carefully  cleaned,  and  the  gage  surface  brought  to  approximately  the 
same  degree  of  reflectivity  as  in  the  previous  run.  The  tube  was  then  evac¬ 
uated  down  to  about  10  -  20  microns  pressure.  The  leak  rate  was  then 
measured.  If  the  rate  did  not  exceed  10  microns  per  minute,  the  tube  was 
considered  to  be  in  satisfactory  condition. 

Next  the  tube  was  filled  with  the  test  gas  to  the  required  test  pressure 
Pj  (normally  either  1  mm  or  5  mm  Hg).  To  further  minimize  the  contami- 
nation  of  the  test  gas  by  impurities  due  to  outgassing  of  the  walls,  the  test 
gas  mixture  was  allowed  to  flow  through  the  tube  at  a  rate  of  approximately 
20  tube  volumes  per  minute  for  about  20  -  30  minutes.  The  initial  pressure 
in  the  driven  tube  was  measured  with  a  McLeod  gage. 

The  driven  section  of  the  2"  shock  tube  is  equipped  with  six  port 
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holes  equally  spaced  at  11"  apart  as  shown  in  Figure  3.  5.  The  model  nor¬ 
mally  is  located  about  1/8  in.  behind  position  6.  The  signal  from  the  piezo¬ 
electric  pressure  gage  located  at  position  1  was  used  to  trigger  the  scopes 
displaying  heat  transfer  gage  output. 

Shock  velocity  was  monitored  by  measuring  the  elapsed  time  for  the 
incident  shockwave  to  travel  between  positions  3-5,  4-5  and  5-6.  Pressure 
gages  were  used  as  time  of  arrival  sensors.  The  last  two  intervals  were 
used  to  determine  the  shock  velocity  for  the  particular  test  run.  The  scatter 
between  data  for  these  two  intervals  is  invariably  within  the  reading  accur¬ 
acy  of  the  wave  transit  time  for  one  interval  (+  1-1/2%).  This  indicates  the 
existence  of  only  small  incident  shock  wave  attenuation  effects  (if  any).  A 

two  color  pyrometer  was  used  to  monitor  the  luminosity  behind  the  incident 

o 

shock  wave.  The  violet  range  extended  between  3500  -  4500  A,  while  the 

o 

red  region  response  was  between  6000  -  8000  A.  The  history  of  the  pressure 
behind  the  incident  shock  wave  was  also  observed  by  displaying  the  output 
of  a  pressure  gage  from  position  5.  The  output  was  filtered  to  eliminate  the 
ringing  of  the  gage  whose  natural  frequency  lies  around  40  kc. 

The  heat  transfer  rate  was  determined  with  a  calorimeter  gage 
mounted  at  the  stagnation  point  of  the  hemispherical  model  1/2  in.  in  diam¬ 
eter  (Figure  3.6).  The  design  dimensions  of  the  gage  were  1/16  in.  wide  by 
3/  16  in.  long.  The  gage  element  itself  was  made  of  0.  004  in.  thick  Hytemco 
(iron-nickel)  alloy.  The  leads  were  brought  out  from  the  gage  through  the 
internal  cavity  of  the  model  and  its  supporting  sting  and  were  connected  in 
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a  circuit  shown  schematically  in  Figure  3.7.  One  of  the  features  (Ref.  17) 
of  the  circuit  is  the  grounding  of  the  instrument  circuit  in  common  with  the 
driver  circuit  and  the  connection  of  the  gage  leads  across  a  balanced  differ¬ 
ential  amplifier.  These  steps  allow  extraneous  electrical  effects  observed 
earlier  in  the  experiments  to  be  reduced  to  a  level  such  that  they  do  not 
preclude  the  observation  of  the  heat  transfer  signal. 

The  output  of  the  heat  transfer  gage  was  displayed  and  photographed 
on  four  oscilloscopes  in  parallel,  each  with  a  different  sweep  rate.  The 
rates  were  20,  10,  5  and  2^sec/cm.  Since  the  region  of  interest  (corre¬ 
sponding  to  the  steady  flow  conditions)  lies  approximately  within  the  first 
20  H  sec  from  the  moment  of  the  bow  shock  formation,  either  the  five  or 
two ^ sec/cm  records  were  used  to  read  off  the  heat  transfer  rate. 

The  relation  between  the  heat  rate  and  the  change  of  the  voltage  with 
time  across  the  gage  is  expressed  by  the  simple  expression 


and 


where 


q 


dT 

dt 


K 


pc  L 
a 


p  =  density  of  gage  material 
C  =  specific  heat  of  gage 
L  =  gage  thickness  ^ 

a  =  temperature  coefficient  of  resistivity 
=  initial  voltage  across  gage 
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Since  the  bulk  properties  of  the  gage  material  are  well  known,  the 
coefficient  K  can  be  easily  calculated  and  the  gage  response  established 
without  an  extensive  calibration  procedure.  The  temperature  coefficient  of 
resistivity,  a,  has  been  checked  experimentally,  however,  in  previous 
studies. 

Stagnation  point  radiant  heat  energy  flux  has  been  measured  by  means 
of  a  newly  devised  technique  which  involves  the  use  of  a  cavity  gage  as  shown 
in  Figure  3.  8.  Depending  on  the  proper  design  of  the  entrance  and  the  loca¬ 
tion  of  the  cavity  gage  within  the  hemispherical  model,  almost  all  of  the 
radiant  energy  originating  from  a  known  volume  of  the  highly  heated  gas  cap 
and  incident  upon  the  cavity  entrance  slit  is  absorbed  within  the  internal 
walls  of  the  gage.  In  addition,  the  electrons  generated  by  high  energy  pho¬ 
tons  will  be  prevented  from  escaping.  The  energy  entering  the  cavity  gage 
will  be  absorbed  by  the  thin  (about  0.  2^)  platinum  film  deposited  on  its 
internal  surfaces.  By  segmenting  the  gage  as  shown  in  Figure  3.  8  the  one- 
dimensionality  of  its  operation  is  preserved  and  its  sensitivity  thereby  in¬ 
creased  considerably.  Due  to  limited  size  of  the  model  for  this  study  (as 
dictated  by  the  blockage  effect  in  the  2"  internal  diameter  shock  tube)  the 
test  configuration  allowed  only  the  measurement  of  radiation  along  the  stag¬ 
nation  streamline  with  a  gage  arrangement  shown  in  Figure  3.  9. 

3.4  Shock  Tube  Parameters 

Using  the  equilibrium  thermodynamic  properties  computed  by 
Browne  (Ref.  5)  for  the  established  two  gas  mixtures,  one  with  the  nominal 
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content  of  9%  CO2,  90%  N2  and  1%  A  and  the  other  with  25%  CO2.  74%  N2 
and  1%  A,  various  gas  dynamic  properties  associated  with  traveling  and 
stationary  shock  fronts  were  calculated  and  charted  as  functions  of  shock 
velocity  and  the  driven  tube  initial  pressures. 

The  tabulated  values  of  the  equilibrium  thermodynamic  properties 
were  plotted  and  the  solution  of  the  Rankine-Hugoniot  shock  relations  ob¬ 
tained.  All  calculations  were  made  for  three  initial  driven  tube  pressures 
of  =  mm,  5  mm  and  25  mm  of  Hg.  The  initial  driven  tube  temperature, 
Tp  was  taken  to  be  294°K.  The  range  of  shock  velocities,  U^,  between 
10,  000  ft/sec  and  30,000  ft/ sec  was  covered.  However,  due  to  the  limited 
range  of  the  available  thermodynamic  properties  which  extended  only  up 
to  15,  000°K  some  of  the  shock  tube  parameters  could  not  be  calculated  for 
all  velocities. 

a.  Incident  Shock 

The  state  of  the  gas  on  both  sides  of  a  normal  shock  is  related  by 
three  basic  equations  corresponding  to  the  physical  principles  of  the  con¬ 
servation  of  mass,  momentum  and  energy.  These  are; 
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where  h 


=  enthalpy 


p  =  pressure 
u  =  velocity 
p  =  density 

and  subscripts  1  and  2  refer  to  regions  ahead  and  behind  the  incident  shock. 
The  simultaneous  solution  of  these  equations  together  with  the  available 
tabulation  of  the  equilibrium  state  produces  results  which  are  plotted  in 
Figures  3.  10,  3.  13,  3.  16,  3.  19  and  3.  21  for  the  9%  CO2  mixture  and  in 
Figures  3.  22,  3.  25,  3.  28,  3.  31  and  3.  33  for  the  25%  CO^  mixture, 
b.  Stationary  Normal  Shock  (Laboratory  Coordinates) 

After  the  passing  of  the  incident  shock  the  model  finds  itself  in  a 
supersonic  quasi-steady  flow  with  the  properties  of  region  2  ahead  of  the 
model  bow  shock  wave.  The  flow  close  to  the  axis  of  symmetry  is  decelerated 
to  subsonic  velocity.  The  state  of  the  gas  along  the  stagnation  streamline  on 
both  sides  of  the  bow  shock  is  governed  by  the  same  three  conservation  laws. 


P2U2  =  P3U3 


where  subscript  3  refers  to  the  region  behind  the  bow  shock.  The  calcula¬ 
tion  produced  the  results  which  are  presented  in  Figures  3.  11,  3,  14,  3.  17, 
3.  20  for  the  9%  CO2  mixture  and  in  Figures  3.  23,  3,  26,  3.  29  and  3,  32  for 


the  25%  CO2  mixture. 
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c.  Stagnation  Region 

The  flow  along  the  stagnation  streamline  in  region  3  is  decelerated 

to  a  stop  at  the  body.  The  stagnation  conditions  are  found  as  follows;  The 

2 

enthalpy  is  given  directly  by  h  =  h^  +  1/2  .  The  stagnation  pressure 

is  computed  approximately  by  assuming  incompressible  flow  + 

2 

1/2  .  Both  the  stagnation  temperature  and  the  corresponding  density 

are  read  off  directly  from  the  thermodynamic  state  charts  and  plotted  in 
Figures  3.  12,  3.  15  and  3.  18  for  9%  CO2  mixture  and  in  Figures  3.  24,  3.  27 
and  3.  30  for  the  25%  CO2  gas. 

d.  Composition  of  Stagnation  Region  Gas 

The  equilibrium  concentration  of  the  various  species  appearing  in 
the  dissociated  and  ionized  gas  when  heated  to  the  enthalpy  level  correspond¬ 
ing  to  the  stagnation  point  had  also  been  obtained  and  plotted  as  a  function 
of  shock  velocity.  These  are  shown  in  Figures  3.  34  through  3,  37  for  both 
the  9%  and  25%  CO2  mixtures  and  for  the  two  test  cases  of  P^  =  1  and  5  mm 
Hg.  Electron  concentration  for  pure  and  pure  are  also  included  in 
Figures  3.  38  and  3.  39  respectively. 

e.  Radiant  Heat  Flux  to  the  Stagnation  Point 

The  radiation  heat  flux  to  the  stagnation  point  of  the  1/2  in.  diameter 
model  was  calculated  using  appropriate  emissivities  obtained  for  the  heated 
atmospheres  of  Venus  (Ref.  18).  The  stagnation  point  region  was  replaced 
by  an  infinite  slab  at  a  uniform  temperature,  T  ,  and  density,  P  ,  whose 

9  8 

thickness,  6 ,  is  the  shock  stand  off  distance  given  by 
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6  ^  2/3 


where  =  radius  of  the  hemispherical  model. 

The  calculations  were  done  for  both  atmospheric  models  and  for 
three  driven  tube  initial  pressures.  The  results  are  presented  in  Figures 
3.  40  and  3.  41. 

3.  5  Results  and  Discussion 

In  the  simulation  of  hypersonic  flight  situations  at  the  stagnation 
point  of  a  model  in  the  shock  tube  the  corresponding  stagnation  enthalpy  and 
pressure  are  produced  in  a  different  manner  than  in  the  flight  case.  The  test 
gas  is  first  compressed  and  accelerated  by  the  incident  shock  wave  generated 
when  the  diaphragm  dividing  the  driver  and  the  driven  tube  is  ruptured.  The 
flow  behind  the  shock  front,  although  moving  with  high  velocity,  is  only 
slightly  supersonic  due  to  the  high  temperature  or  high  sonic  velocity  in  the 
gas  processed  by  the  advancing  incident  shock  wave.  This  gas  will  now  flow 
around  our  model.  Since  a  blunt  body  is  being  investigated,  a  detached  bow 
wave  will  be  established.  The  gas  along  the  stagnation  streamline  will  under- 
go  compression  essentially  across  a  normal  shock  and  finally  will  be  stag¬ 
nated  at  the  surface  of  the  body. 

As  long  as  the  boundary  layer  thickness  is  small  in  comparison  with 
the  shock  layer  whose  measure  at  the  stagnation  point  is  the  detachment  dis¬ 
tance,  the  boundary  layer  problem  can  be  considered  separately  from  the 
inviscid  flow  field  and  hence  independently  of  the  Reynolds  number.  The 
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only  requirements  for  accurate  simulation  are  the  stagnation  enthalpy,  the 
velocity  gradient  at  the  stagnation  point,  and  gas  equilibrium  in  the  various 
flow  processes.  For  the  hemispherical  body  the  Newtonian  approximation 
predicts  that  the  velocity  gradient  is  independent  of  Mach  number.  This  has 
been  confirmed  experimentally.  The  correlation  for  the  heat  transfer  rate 
through  the  laminar  boundary  layer  at  the  stagnation  point  can  be  written 


where  is  body  nose  radius 

h^  is  wall  enthalpy 

Thus,  the  environment  created  in  a  shock  tube  can  be  directly  applied  to  the 
solution  of  the  stagnation  point  heating  problem. 

The  performance  of  the  shock  tube  in  the  mixtures  of  CO^.  N2  and 
A  does  not  differ  appreciably  from  its  performance  in  air.  Maximum  shock 
velocities  obtained  during  this  experimental  program  for  values  of  5  mm 
and  1  mm  Hg  were  approximately  24,  500  ft/ sec  and  30,  200  ft/ sec.  In  both 
cases  maximum  capacitor  energy  <about  50,  000  joules)  was  discharged  in 
the  driver  for  heating  the  helium. 

The  constituent  concentrations  of  the  mixtures  tested  differed  slightly 
from  the  nominal  values  listed  earlier. 

Nominal  Actual 

1.  9%C02,  90%  N2,  1%A  1.  8.80%CX)2,  89.98%N2,  1.22%  A 

2.  25%  CO2,  74%  N2,  1%  A  2.  25.  10%  CO2,  73.  78%  N2,  1.  12%  A 
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The  test  gas  mixtures  were  considered  to  be  sufficiently  close  to  the  specified 
values  to  satisfy  the  objectives  of  the  study. 

During  each  run,  wall  pressures  histories  behind  the  incident  shock 
wave  were  recorded.  A  typical  oscillogram  of  signal  from  a  pressure  gage 
located  at  position  five  is  shown  in  Figure  3.42,  corresponding  to  a  shock 
velocity  =  20,  300  ft/  sec  and  initial  tube  pressure  of  5  mm  Hg.  The  effect 
of  boundary  layer  growth  behind  the  incident  wave  is  indicated  by  the  side 
wall  pressure  rise  with  time.  Note  that  this  rise  is  small  during  the  passage 
of  the  test  gas,  however.  Pressure  forces  of  the  type  shown  in  Figure  3.42 
are  complicated  by  large  oscillations  caused  by  the  response  of  the  gage  and 
associated  electronic  circuits.  They  are  used  primarily  to  verify  the  test 
pressure  level  and  to  qualitatively  interpret  the  phenomena  occurring  in  the 
tube  process  (arrival  of  expansion  wave  •  not  seen  in  Figure  3.  42  -  and  boun¬ 
dary  layer  effect  on  pressure  level,  for  example). 

The  heat  transfer  measurements  were  made  at  the  stagnation  point 
of  a  0.  5  in.  diameter  hemisphere -cylinder  model  using  a  calorimeter  heat 
transfer  gage  (Ref.  19).  The  gage  material  used  for  the  0.5  in.  diameter 
model  was  .  004  in.  thick  Hytemco.  A  few  experiments  were  performed 
with  a  0.  75  in.  diameter  model  and  the  gage  used  0.  003  in.  thick  platinum. 
According  to  the  analytical  results  of  Ref.  19,  these  thicknesses  can  only 
introduce  the  order  of  1%  error  in  the  measured  heating  results.  Their 
use  is,  therefore,  justified  by  practical  considerations.  The  oscillogram 
of  the  signals  from  both  type  gages  are  shown  in  Figures  3.44a  and  3.44b 
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respectively.  The  conditions  of  both  tests  were  very  similar.  The  shock 
velocities  were  29,200  and  29,600  ft/ sec,  and  the  initial  pressures  1  mm 
Hg  in  both  cases,  '^he  first  test  was  in  25%  CO^  mixture;  the  second,  in 
air.  Although  in  this  particular  case  the  noise  level  of  the  response  of 
either  gage  seems  to  be  of  the  same  magnitude,  closer  examination  of  the 
traces  w'ill  prove  that  the  platinum  gage  response  is  cleaner.  This  was 
also  observed  in  other  runs.  The  initial  jump  of  the  DC  level  with  the 
Hytemco  gage  was  also  observed  to  occur  with  the  platinum  gage.  Lack  of 
similarity  of  the  precursor  response  during  the  two  or  three  microseconds 
before  shock  arrival  cannot  be  explained  at  present. 

The  oscillogram  from  a  calorimeter  heat  gage  at  the  stagnation 
point  of  the  0.  5  in,  model  is  shown  in  Figure  3.  43.  The  instantaneous  rate 
of  change  of  the  output  together  with  the  known  properties  of  the  Hytemco 
alloy  gives  the  unknown  heat  rate  for  a  given  time  interval.  A  corresponding 
trace  of  the  optical  pyrometer  sensing  the  radiation  across  the  incident 
shock  of  station  5  is  also  shown.  Taking  the  time  interval  during  which  the 
intensity  of  radiation  is  level  and  allowing  for  the  formation  of  the  bow  shock 
(2  to  3 /is),  the  useful  test  time  was  arrived  at  and  the  slope  of  the  gage 
signal  during  that  period  was  measured,  giving  the  unknown  heat  transfer 
rate. 

The  experimental  study  involved  a  series  of  runs  in  the  9%  and  25% 
CO2  mixtures  over  a  wide  shock  velocity  range  (note  that  does  not  equal 
U£),  simulating  hypersonic  entry  at  altitudes  between  240,  000  and  260,  000 
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ft.  of  the  atmospheric  model  I.  The  results  are  shown  in  Figure  3.45. 

Points  marked  by  circles  and  triangles  are  from  tests  in  9%  and  25%  CO^ 
mixtures,  respectively.  In  addition,  some  data  obtained  in  pure  CO2  and 
N2  are  included.  The  limited  results  in  CO^  heat  transfer  measurements 
of  Rutowski  and  Chan  (Ref.  20)  at  low  enthalpies  are  also  entered. 

The  theoretical  curves  representing  Scala'  s  calculation  of  equilibrium 
boundary  layer  for  dissociated  and  ionized  nitrogen  are  drawn  together  with 
his  early  estimated  region  of  uncertainty  for  partially  ionized  carbon  dioxide. 

At  (h  -  h  )  values  below  24,000  BTU/lb  the  experimental  data  follow  the 

s 

theoretical  curve.  At  higher  enthalpies  the  present  data  give  heat  transfer 
rates  somewhat  lower  than  the  ionized  theory  for  N2.  The  pyrometer  records 
indicate  that  the  simulation  of  flight  velocities  below  35,  000  ft/ sec  gives 
valid  data  with  sufficient  test  time.  The  results  corresponding  to  higher 
flight  velocities  may  have  insufficient  test  time  for  accurate  readings.  This 
limitation  is  more  likely  to  exist  for  the  CO2  mixtures  than  for  air  and  nitro¬ 
gen  because  of  the  apparent  long  relaxation  times  for  the  non-equilibrium 
process  behind  the  incident  shockwave  (see  Figure  3.43).  Such  an  effect 
could  cause  the  observation  of  values  lower  than  actually  exist  and  is.  there¬ 
fore,  Quite  important.  This  possible  restriction  on  the  results  is  imposed 
by  the  present  2"  test  facility  and  will  be  remedied  by  the  installation  of  the 
6"  shock  tube  scheduled  for  operation  shortly. 

It  is  interesting  here  to  note  that  the  electron  densities  of  the  stagna¬ 
tion  point  as  calculated  for  9%  and  25%  CO2  mixtures,  pure  N2  and  pure  CO2 
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differ  little  from  each  other;  this  can  be  seen  in  Figures  3.  34  through  3.  37. 

The  contribution  of  radiation  to  the  measured  heat  transfer  rate  was 
calculated  and,  in  the  case  of  both  gas  mixtures,  its  magnitude  is  less  than 
8%  assuming  complete  absorption  by  the  gage  surface.  This  is  the  theoreti¬ 
cal  upper  limit  of  possible  radiative  error  in  the  measured  boundary  layer 
heating.  Figure  3.  46  shows  the  response  of  the  cavity  radiation  gage  during 
a  test  run  in  the  25%  CO2  mixture.  It  is  pointed  out  that  the  sensing  ele¬ 
ment  in  the  cavity  gage  is  a  thin  film  surface  temperature  resistance  ther¬ 
mometer  so  that  its  responseis  different  from  that  of  a  calorimeter  gage. 
For  example,  for  a  constant  heating  pulse  applied  in  a  step  function  at  t  s  0, 
the  thin  film  response  will  be  proportional  to  VT while  the  calorimeter 
response  will  be  proportioned  to  t.  The  signal  shown  exhibits  some  elec¬ 
tronic  noise  just  prior  to  shock  arrival  over  a  time  period  similar  to  that 
observed  with  the  calorimeter  gage  traces.  Again,  this  is  undoubtedly  due 
to  rapid  charge  equilibration  about  the  model  in  the  ionized  flow.  The  physi¬ 
cal  mechanisms  causing  these  early  signals  are  not  clear;  however,  it  again 
appears  that  they  become  unimportant  soon  after  the  shock  wave  reaches 
the  model.  The  relatively  clean  response  of  the  gage  during  approximately 
the  lOp  sec  following  the  shock  wave  arrival  indicates  that  this  is  the  time 
for  passage  of  the  shocked  test  gas.  The  subsequent  fall  and  rise  are  re¬ 
sponses  to  the  irregular  driver  gas  flow  that  then  impinges  on  the  model. 
Note  that  a  drop  in  response  is  an  indication  of  suddenly  reduced  heating 
rate.  The  radiative  heating  rate  at  the  stagnation  point  corresponding  to 
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the  gage  signal  shown  in  Figure  3.  46  has  not  yet  been  reduced  accurately. 
However,  an  approximate  reduction  shows  it  to  be  within  30<%  of  the  calcu¬ 
lated  value.  Therefore,  the  conclusion  is  drawn  that  the  radiation  contribu¬ 
tion  to  measured  heating  rate  values  is  small  (10%  or  less).  This  conclusion 
is  now  based  upon  both  experimental  and  theoretical  grounds. 

3.  6  Conclusions 

The  experimental  program  of  stagnation  point  heat  transfer  rates 
measurement  on  a  blunt  body  in  two  simulated  atmospheres  of  Venus  supplied 
data  which  indicate  that  within  the  experimental  accuracy  no  discernible 
difference  between  the  heat  transfer  rate  in  9%  CO2  mixture  and  in  the  25% 
CO2  mixture  can  be  observed.  For  conditions  corresponding  to  flight  velo¬ 
cities  below  35,  000  ft/  sec  the  experimental  data  agree  well  with  theory. 

For  higher  values  of  (h^  -  h^)  or  corresponding  flight  velocities  the  test 
results  may  suffer  due  to  short  test  time  caused  by  slow  approach  to  equili¬ 
brium  of  reaction  processes  behind  the  normal  shock.  The  new  6"  shock 
tube  now  in  process  of  installation  will  increase  considerably  the  available 
test  time.  Calculation  of  the  upper  limit  of  equilibrium  radiation  from  the 
shock  layer  to  the  stagnation  point  indicates  that  its  contribution  to  the  mea¬ 
sured  rates  is  small.  Measurement  of  the  radiation  by  the  newly  developed 
total  radiation  gage  gives  values  close  to  the  theoretical  prediction. 

Further  studies  of  equilibrium  and  non-equilibrium  radiation  and  the 
influence  of  non-equilibrium  flows  on  the  convective  heat  transfer  rate  are 
recommended.  The  behavior  of  the  laminar  and  turbulent  boundary  layers 
and  the  transition  phenomena  should  also  be  studied  in  the  hypervelocity  shock 
tube  facility. 
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Figure  3.  1  Design  of  Arc  Heated  Shock  Tube. 
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Figure  3.  2  Typical  Electrical  Characteristics  of  Prototype  Shock  Tube 
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Figure  3,  3  High  Enthalpy  Shock  Tube  Facility 
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Figure  3.4  Arc  Heated  Helium  Driver  Shock  Tube,  Mg  vs. 

(2"  Shock  Tube  -  Max.  Available  Energy;  47,000  Joules) 
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Figure  3.  6  Half -inch  Model  with  Heat  Transfer  Gage 
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b.  SCHEMATIC  OF  GAUGE  CIRCUIT 


Figure  3.  7  Heat  Transfer  Rate  Measurement  Set-up. 
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ADVANTAGES: 

1.  WINDOWLESS  AND  NO  SURFACE  EFFECT  UNCERTAINTY.  ACCEPTS  AND 
MEASURES  ALL  INCIDENT  RADIATION 

2.  COLUMATION  ALLOWS  MEASUREMENT  OF  TOTAL  EMISSION  FROM  LOCAL 
AREAS  OF  HEATED  GAS  WITH  Us  RESOLUTION. 


SCHEM.  GAGE  CIRCUIT 


MOVING  SHOCK  WAVE 


Figure  3.  8  Cavity  Gage  for  Total  Radiative  Heat  Transfer  Measurements 
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Figure  3.  9 


Figure  3.  10  Temperature  Behind  Incident  Shock 
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Figure  3.  11  Temperature  Behind  Normal  Shock 


Figure  3.  12  Stagnation  Temperature 
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Figure  3.  14  Pressure  Behind  Normal  Shock 


Figure  3.  15  Stagnation  Pressure 


Figure  3.  16  Enthalpy  Behind  Incident  Shock 


Figure  3.  17  Enthalpy  Behind  Normal  Shock 


Figure  3.  18  Stagnation  Enthalpy 
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Figure  3.  19  Density  Behind  Incident  Shock 


Figure  3.  20  Density  Behind  Normal  Shock 
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Figure  3,21  Velocity  Behind  Incident  Shock 
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Figure  3,  22  Temperature  Behind  Incident  Shock 
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Figure  3.  23  Temperature  Behind  Normal  Shock 
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Figure  3.  24  Stagnation  Temperature 
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Figure  3.  27  Stagnation  Pressure 
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Figure  3,  28  Enthalpy  Behind  Incident  Shock 
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Figure  3.  29  Enthalpy  Behind  Normal  Shock 
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Figure  3.  30  Stagnation  Enthalpy 
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Figure  3,  31  Density  Behind  Incident  Shock 
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Figure  3.  32  Density  Behind  Normal  Shock 
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Figure  3.  33  Velocity  Behind  Incident  Shock 
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Figure  3.  34  Equilibrium  Composition  of  Venus  Atmosphere  at  Stagnation 
Point  of  a  Hemispherical  Model 
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Figure  3.  35  Equilibrium  Composition  of  Venus  Atmosphere  at  Stagnation 
Point  of  a  Hemispherical  Model 
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Figure  3.  36  Equilibrium  Composition  of  Venus  Atmosphere  at  Stagnation 
Point  of  a  Hemispherical  Model 
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Figure  3.  37  Equilibrium  Composition  of  Venus  Atmosphere  at  Stagnation 
Point  of  a  Hemispherical  Model 
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Figure  3,  38  Equilibrium  Concentration  of  Electrons  at  the  Stagnation 
Point  of  a  Hemispherical  Model 
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Figure  3.  39  Equilibrium  Concentration  of  Electrons  at  the  Stagnation 
Point  of  a  Hemispherical  Model 
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Figure  3.40  Equilibrium  Radiative  Heat  Transfer  to  the  Stagnation  Point 
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Figure  3.  42  Oscillogram  of  Pressure  Behind  Incident  Shock 
(Signal  Unfiltered) 
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Figure  3.  43  Typical  Oacillograma  from  Heat  Transfer  Gage  and  Corre> 
spending  Pyrometer  Signal  of  Incident  Shock  Wave 
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Figure  3.  44a  Typical  Oscillogram  from  Hytemco  Calorimeter 
Heat  Transfer  Gage 
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Figure  3.  44b  Typical  Oscillogram  from  Platinum  Calorimeter  Heat 
Transfer  Gage 
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Arrival  of  Shock  Wave 
Test  Time 


Figure  3.  46  Total  Radiation  Thin  Film  Cavity  Gage  Response 


III-4  LARGE.  HIGH  ENTHALPY  SHOCK  TUBE  DEVELOPMENT* 


Using  as  the  basis  the  2"dia.  arc -heated  shock  tube,  a  larger  facility 
capable  of  producing  shock  velocities  up  to  40,  000  ft/ sec  was  designed.  The 
general  arrangement  is  shown  in  Fig.  4.  1.  The  driver  was  made  larger  by 
increasing  its  length  to  4  1/2  ft.  and  its  diameter  to  3  1/2  in.  The  energy 
stored  in  capacitors  used  to  heat  the  helium  was  also  increased  to  320K 
joules  by  adding  several  more  capacitor  banks.  Most  of  the  basic  design 
features  of  the  prototype  unit  were  retained.  The  arc  is  initiated  by  a  tung¬ 
sten  wire  suspended  coaxially  inside  a  Lexan  insulator  between  the  main 
electrode  and  the  spider  at  the  diaphragm  end  of  the  driver.  This  assures 
uniform  wtdl  losses  along  the  driver.  To  keep  the  energy  dissipation  time 
short  (about  50  microseconds)  a  parallel- series  arrangement  of  the  capaci¬ 
tor  circuit  is  employed.  The  driver  tube  is  30  ft.  long  and  6  in.  inside 
diameter.  It  is  made  of  steel  permitting  the  mounting  of  pressure  gages 
along  the  tube  to  evaluate  flow  properties  and  to  sense  the  arrival  of  driver 
expansion  wave.  Other  instrumentation  such  as  the  cavity  gage  which 
requires  contact  with  the  test  gas  can  also  be  used.  The  tube  is  relatively 
easy  to  clean  and  is  capable  of  withstanding  the  high  pressure  experienced 
at  high  P^  values. 

Considerable  improvement  of  the  range  of  test  conditions  is  achieved 
<Principal  Investigators:  J.  S.  Gruszczynski,  C.  J.  Harris,  D.  A.  Rogers 
and  W,  R.  Warren. 
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by  the  larger  tube  in  comparison  to  the  2  in.  prototype  facility.  This  is 
mainly  reflected  by  the  test  time  available  at  a  given  initial  driven  tube 
pressure  and  shock  velocity.  Figures  4.  2  and  4.  3  illustrate  this  point. 
The  range  of  test  conditions  behind  incident  shock  wave  as  shown  in  Fig. 

4.2  is  limited  by  the  maximum  energy  input  into  the  driver  corresponding 
to  the  final  pressure  in  the  driver  of  2000  atmospheres  and  8000°K  tem¬ 
perature  on  one  side,  and  by  the  test  time  which  precludes  the  operation 
below  100  microns  on  the  other  side.  A  typical  low  W/Cj^A,  low  tra¬ 
jectory  is  superimposed  upon  the  shock  tube  performance  together  with 
the  altitude  scale  of  a  model  of  Venusian  atmosphere. 

Fig.  4.  3  shows  the  performance  of  the  shock  tube  related  to  the  sim¬ 
ulation  of  stagnation  point  conditions.  Here  we  see  that  the  flight  velocities 
of  40,  000  ft/ sec  can  easily  be  simulated  for  stagnation  pressures  up  to  100 
atm. 

Testing  has  not  commenced  in  this  facility  as  yet.  All  the  hardware 
has  been  procured  and  is  now  installed.  Fig.  4.  4  shows  the  partially  com¬ 
pleted  capacitors  room  with  the  driver  in  place.  The  driven  tube  is  shown 
in  Fig.  4.  5.  The  test  section  which  will  be  mounted  on  the  end  of  the  driver 
tube  is  not  shown. 

The  full  facility  has  begun  initial  operation  and  is  being  checked  out. 
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Figure  4.  2  Performance  Curve  for  Hypervelocity  Shock  Tube 


Figure  4.  3  Performance  Curve  for  Hypervelocity  Shock  Tube 


Fiaure  4.  ^  Driven  Tube  Hypervelocity  Shock 


III-5  HEAT  PROTECTION  STUDIES* 


5.  1  Introduction 

The  thermal  protection  by  ablation  of  vehicles  entering  simulated 
planetary  atmospheres  was  the  subject  of  the  present  study.  The  capability 
for  such  a  study  was  provided  by  the  Space  Sciences  Laboratory' s  tandem 
Gerdien  arc  --a  high  purity,  high  enthalpy  design.  The  tandem  Gerdien 
design  incorporating  the  latest  "state  of  the  arc"  features  was  adapted  to 
free  jet  testing.  Well  known  ablation  materials  typical  of  three  classes  of 
ablators  were  tested  in  the  arc-heated  gas  mixtures  representative  of  Mar¬ 
tian  and  Venusian  atmospheres.  The  performance  of  each  material  in  each 
gas  mixture  is  evaluated  and  compared  with  similar  performance  in  air. 

The  results  illustrate  the  severity  of  the  entry  problem  for  Mars  and  Venus 
relative  to  Earth  re-entry  for  high  heating  rate  or  primarily  steady  state 
ablation  situations.  The  effects  of  non-steady  ablation  and  gas  radiation  on 
material  performance  were  not  studied. 

5.  2  Arc  Test  Facility 

The  tandem  Gerdien  arc  design  was  adapted  to  free  jet  testing  of  abla¬ 
tion  materials.  The  arc  heater  provides  contaminant -free  high  enthalpy 
flow  due  to  its  constrictor  plate  design  and  the  flow  velocity  distribution  in 
the  arc  (Figure  5.  1).  Argon  sheathed  electrodes  were  incorporated  in  order 
to  substantially  reduce  the  electrode  consumption  rate.  Since  the  inlet  gas 
is  divided  into  four  portions,  two  flowing  through  the  arc,  plenum  chamber 
and  nozzle  (test  gas)  and  two  flowing  past  the  electrodes  to  a  dump  valve, 
"iPrincipal  Investigators:  N.  S.  Diaconis,  R.  A.  Sheridan  and  W.  R.  Warren 
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no  electrode  materiale  or  argon  ia  entrained  in  the  test  gas. 

The  tandem  Gerdien  arc  was  developed  by  McGinn  at  the  Space  Sciences 
Laboratory.  It  incorporated  the  principles  and  experience  gained  in  the 
seven  years  of  arc  technology  at  the  Laboratory.  The  arc  unit  included  the 
latest  feature  of  arc  development  --  argon  swirl  electrodes.  The  graphite 
electrodes  are  shrouded  in  argon  contained  by  a  water-cooled  copper  housing. 
The  argon  prevents  combustion  of  the  electrodes,  and  the  argon  and  carbon 
are  dumped  to  prevent  contamination  of  the  test  gas. 

Other  modifications  of  the  arc  unit  were  needed  to  adapt  it  to  free  jet 
testing.  The  model  actuation  mechanism  is  shown  in  Figure  5.  2.  A  slot 
was  made  in  the  plenum  wall  for  full  viewing  of  the  models;  the  back  pots 
were  modified  to  accommodate  the  argon  swirls,  and  other  changes  were 
made  to  facilitate  the  operation  of  the  arc.  Power,  water,  gas  and  hydrau¬ 
lic  lines  were  installed  along  with  sufficient  instrumentation  to  monitor  the 
operational  parameters  of  the  arc. 

An  advantage  gained  when  the  arc  became  operational  was  the  utiliza¬ 
tion  of  the  new  power  source.  The  rectified  power  supply  provides  a  more 
constant  source  than  the  old  supply  and  at  higher  power  levels.  Improved 
operation  of  the  arc  was  clearly  noted. 

After  the  initial  trial  runs  with  the  newly  built  free  jet  arc  (Fig.  5.  3) 
the  purity  of  the  flow  was  checked  spectrographically.  The  amount  of  inlet 
flow  past  the  argon  sheathed  electrodes  and  the  amount  exhausted  were 
varied  until  the  test  gas  species  showed  a  high  degree  of  purity  (contamina- 
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tion  approximately  0.  1%  by  mass). 

Some  difficulties  were  encountered  in  the  initial  operation  of  the  arc 
unit.  Water  leaks  and  arc  covers  to  the  arc  parts  necessitated  several  de¬ 
lays  while  repairs  were  made.  The  arc  was  debugged  and  brought  under 
control  and  the  performance  was  measured  to  a  satisfactory  degree.  The 
arc  current  and  voltage,  the  plenum  pressure,  the  purity  of  the  flow,  the 
stagnation  temperature  and  stagnation  enthalpy  were  determined  for  the  arc. 

The  stagnation  enthalpy  of  the  gas  mixtures  was  determined  with  a 
water  cooled  total  calorimeter  (Figure  5.  4).  The  calorimeter  collects  the 
gas,  cools  it  and  then  exhausts  it  through  a  flow  meter.  The  arc  conditions 
are  maintained  constant  until  the  calorimeter  comes  to  equilibrium.  The 
enthalpy  is  determined  from  the  water  temperature  rise  and  the  weight  flow 
of  gas. 

Model  heat  transfer  rates  are  determined  through  transient  backwall 
temperature  measurements  in  a  cylindrical  slug  located  at  the  stagnation 
point  of  a  plastic  model. 

5.  3  Test  Program 

The  objective  of  the  program  was  to  determine  experimentally  the 
behavior  of  typical  thermal  protection  materials  that  might  be  used  for  Mars 
and  Venus  entry  vehicles.  An  approach  that  would  allow  compensation  for 
any  mis -estimate  of  the  percentages  of  the  nitrogen-carbon  dioxide  mix¬ 
ture  for  the  planets  was  used.  The  percentage  of  CO2  (in  a  C02-N2«  mix¬ 
ture)  was  varied  from  0  to  25%.  The  three  mixtures  tested  then  allow  re- 
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suits  to  be  approximated  for  any  mixture  by  interpolation.  Test  results  in 
the  simulated  atmospheres  were  also  compared  with  air  data. 

Nominally,  the  gas  mixtures  selected  for  the  program  were  air,  nitro¬ 
gen,  and  the  "best"  estimate  and  the  extreme  estimate  for  the  Venus  atmo¬ 
sphere.  The  actual  N2-CO2  desired  and  those  obtained  from  the 

vendor  were  as  follows; 

Mixtures  Specified  Mixtures  Supplied 

90%  N2,  9%  CO2.  1%A  88.6%N2,  10%  CO2.  1.4%  A 

74%N2.25%  CO2,  1%A  74.  1%N2,  24.  7%  CO2,  1.2%  A 

These  were  considered  to  be  satisfactory  mixtures  for  the  purposes  of  the 
program. 

The  materials  selected  for  study  are  representative  of  three  some¬ 
what  general  classes  of  ablators; 

Phenolic  Nylon  --  pyrolizes  and  chars 
ATJ  Graphite  --  combusts  at  the  surface 
Teflon  --  sublimes 

To  permit  comparison  of  the  model  tests  in  the  different  gas  mixtures 
it  is  necessary  to  fix  certain  test  parameters.  The  stagnation  pressure 
was  kept  essentially  constant  and  the  stagnation  enthalpy  relatively  constant 
for  each  series  of  test  gases.  In  this  way  the  results  are  comparable  among 
themselves  and  are  representative  of  stagnation  region  ablation  material 
performance  at  a  specific  flight  velocity  and  flight  altitude. 

The  majority  of  the  tests  were  performed  in  the  stagnation  enthalpy 
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range  of  12,  000  to  14,000  BTU/lb.  A  few  tests  were  made  in  a  10,  000  BTU/ 
lb  enthalpy  range.  The  maximiim  simulated  flight  velocity  and  altitude 
values  were  approximately  26,  500  ft/ sec  and  83  KM  (for  a  simulated  Venu¬ 
sian  atmosphere). 

5.  4  Test  Results  and  Discussion 

The  stagnation  enthalpy  values  determined  with  the  total  calorimeter 
are  shown  in  Table  5.  1.  The  tests  were  repeated  at  least  twice  for  each 
gas  mixture  and  the  results  indicated  that  the  arc  flow  conditions  are  highly 
repeatable  (+2%).  Table  5.  1  shows  the  operating  pressure  in  the  plenum, 
the  enthalpy  increase  of  the  test  gas  due  to  arc  heating,  the  stagnation  en¬ 
thalpy  and  the  wall  enthalpies  at  wall  temperatures  of  1500°R  and  5400°R 
for  each  gas  mixture.  Also  shown  are  the  heat  transfer  rates  for  each  gas 
at  the  two  wall  temperature  conditions  of  interest.  The  calculated  equili¬ 
brium  thermochemical  properties  given  in  References  24  and  25  were  used  to 
determine  the  necessary  properties  used  below  in  the  data  presentation 
(hoo,  h^). 

Model  stagnation  point  heat  transfer  rate  data  (non-ablating)  was  ob¬ 
tained  for  each  gas.  The  results  indicated  that  many  additional  runs  would 
have  to  be  made  to  obtain  reliable  mean  heat  transfer  rates  for  each  gas. 

The  reasons  for  this  are  that  at  the  high  test  heat  transfer  rates  (of  the 
order  of  2500  BTU/ft  -sec),  first,  the  plastic  calorimeter  holder  ablates 
quickly  near  the  copper  insert,  and  second,  the  surface  temperature  of  the 
copper  quickly  approaches  its  melting  level.  Thus,  during  the  time  for 
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which  the  calorimeter  responds  accurately  it  is  difficult  to  measure  the 
temperature  slope  data  (with  present  equipment).  The  result  is  an  appre¬ 
ciable  scatter  in  the  reduced  heat  transfer  rates  (+  10  to  15%), 

There  was  insufficient  time  to  conduct  the  number  of  tests  required 
for  each  gas  in  this  program;  therefore,  a  semi -empirical  method  was  used. 
The  mean  heat  transfer  rate  in  air  was  determined  from  several  experiments 
at  the  high  h^  levels  (Table  5.  1).  Theoretical  relationships  for  air,  nitro¬ 
gen,  and  CO2  (Refs.  21  and  22)  were  then  used  to  establish  heat  transfer 
ratios  between  each  test  gas  and  air  at  the  appropriate  stagnation  enthalpy 
level.  Thus,  the  heating  rate  for  each  gas  was  calculated  according  to 
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Since  the  enthalpies  and  theoretical  heating  rates  for  each  gas  differ  from 
the  air  values  only  by  the  order  of  20%  and  10%,  respectively,  this  proce¬ 
dure  should  be  satisfactory.  A  further  justification  is  found  in  the  good 
facility  operation  repeatability.  As  a  check  on  the  measured  heating  rate 
levels,  calculations  for  air  at  the  test  conditions  were  made  for  stagnation 
point  velocity  gradients  corresponding  to  incompressible  potential  flow  and 
the  modified  Newtonian  prediction.  The  potential  value  agrees  (within  3%) 
with  the  mean  measured  value  as  would  be  expected  for  the  subsonic  free 
jet  test  condition  used. 


99 


The  different  ablators  were  tested  for  durations  up  to  6  seconds. 
Longer  test  times  would  be  meaningless  because  of  the  large  amount  of  ma¬ 
terial  loss  and  shape  change  of  the  models.  However,  at  the  heating  rates 
of  interest,  steady  state  ablation  conditions  are  established  for  most  mate¬ 
rials  in  a  small  fraction  of  this  time.  Figure  5.  5  shows  typical  models 
tested  for  four  seconds  duration.  The  shape  changes  are  to  be  noted,  espe¬ 
cially  for  the  teflon.  The  phenolic  nylon  has  a  strong  thick  char  layer  so 
that  the  shape  change  is  not  as  pronounced  as  with  teflon.  However,  the 
shape  change  of  the  virgin  phenolic  nylon  is  quite  large  because  of  the  side 
heating  and  char  growth.  ATJ  graphite *data  is  not  shown,  although  the 
material  was  tested  in  each  gas.  The  mass  loss  in  air  at  a  measured  sur¬ 
face  temperature  of  5400°R  was  significantly  higher  than  expected.  The 
theory  for  graphite  and  its  experimental  verification  are  presented  in  Ref¬ 
erences  Z3  and  24.  The  increased  mass  loss  occurs  after  the  test  flow  is 
terminated,  that  is,  the  surface  continues  to  combust  in  the  room  air  as 
the  model  cools.  Due  to  the  heat  capacity  of  graphite  the  duration  of  tem¬ 
perature  decay  is  several  times  the  test  duration.  This  problem  was  not 
found  for  the  other  materials  tested  since  with  them  little  total  heat  is  ab¬ 
sorbed  during  test,  and  what  is,  stays  near  the  surface  resulting  in  a  short 
cool  down  period.  Thus,  the  graphite  data  are  not  included  since  the  results 
could  be  misleading.  Even  with  the  apparent  (and  erroneous)  increase  in 
mass  loss,  the  material  ablation  rate  was  still  much  less  than  for  PN  at 
the  same  test  conditions.  It  must  be  remembered,  however,  that  graphite 
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blocks  much  less  heat  than  does,  for  example,  PN,  resulting  in  its  high 
internal  temperatures.  Therefore,  mass  loss  becomes  a  relatively  less 
important  design  consideration  for  graphite. 

It  is  noted  that  this  problem  for  graphite  can  be  eliminated  in  future 
studies  by  providing  an  inert  gas  bath  for  the  test  models  during  the  cooling 
period. 

The  models  before  test  (Figure  5.  5)  are  hemisphere-cylinders  1/2  inch 
in  diameter.  A  3/  16  inch  diameter  by  1  inch  long  slug  of  the  same  material 
is  inserted  at  the  stagnation  point  of  the  model.  All  test  measurements  are 
made  for  this  slug.  This  technique  allows  the  obtaining  of  stagnation  region 
data  as  opposed  to  data  integrated  over  the  complete  model  surface.  The 
right  fit  eliminates  effects  due  to  the  model-slug  joint. 

Weight  and  length  measurements  are  taken  from  the  slug  before  and 
after  test.  Typical  test  specimens  are  shown  in  Figure  5.6.  Note  the  char 
cap  positioned  above  the  phenolic  nylon  specimen.  This  becomes  detached 
from  the  slug  when  it  is  removed  from  the  model  after  test.  The  weight 
loss  and  char  thickness  for  the  series  of  runs  for  phenolic  nylon  and  teflon 
are  given  in  Table  II.  The  ratios  of  the  initial  model  radius  (R^  =  1/4") 
to  the  average  radius  after  test  (R^)  are  also  shown. 

Spectrograms  were  taken  of  each  model  in  each  test  gas  to  attempt  to 
justify  the  use  of  a  two-color  pyrometer  to  determine  the  surface  tempera¬ 
ture  history  of  the  test  specimens.  The  spectrograms  showed  that  neither 
color  nor  brightness  temperatures  could  be  determined  with  the  pyrometer 
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for  the  majority  of  the  tests.  That  is,  the  gas  radiation  continuum  was  too 
intense  and  extensive  to  provide  "windows"  of  sufficient  frequency  interval 
to  see  the  model  surface  satisfactorily  with  the  pyrometer.  To  obtain  some 
information  on  the  surface  temperatures,  a  few  additional  spectrographic 
measurements  were  made.  The  temperatures  determined  from  analysis  of 
the  spectrograms  were  as  follows: 

High  Enthalpy  Test  Conditions 

Test  Gas  Material  Surface  Temperature 

10%  CO2,  88.6%  N2,  1.4%  A  AT  J  graphite  3900°K  +  100°K 

Air  ATJ  graphite  2900°K  +  100°K 

Air  Phenolic  nylon  3000°K  +  100®K 

The  surface  temperatures  assumed  for  all  the  tests  were  3000°K  for 
phenolic  nylon  and  1500°R  for  teflon  (based  on  previous  experimental  and 
theoretical  experience  with  teflon  (Refs.  25  and  26).  It  is  interesting  to 
note  that  the  surface  temperatures  of  graphite  are  significantly  higher  in 
the  10%  CO^  mixture  than  in  air.  This  is  expected  since  the  ablation  mech¬ 
anism  for  graphite  is  primarily  combustion.  (Note  that  sublimation  should 
o 

start  at  about  3700  K  for  graphite  at  one  atmosphere. )  It  is  not  anticipated 
that  PN  or  teflon  will  exhibit  this  type  of  change  since  their  ablation  mech¬ 
anisms  are  not  dominated  by  combustion.  However,  this  point  should  be 
checked  in  future  studies. 

Figure  5,  7  shows  the  weight  loss  vs.  time  for  teflon  and  phenolic 
nylon  in  each  test  gas  mixture  at  the  two  nominal  enthalpy  levels.  The  data 
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do  not  show  a  linear  correlation  since  the  model  heating  rate  is  decreasing 
with  time  (nose  radius  changing).  The  data  was  reduced  on  the  basis  of  in¬ 
stantaneous  values.  The  mass  loss  is  computed  from  the  slope  of  the  weight 
loss -time  curve  and  a  heat  transfer  rate,  appropriate  to  the  model  nose 
radius  at  the  same  time  of  test,  is  calculated  according  to 


The  heat  of  ablation  is  defined  as 


M 


♦  *  • 

where  and  M  are  determined  as  above  with  corrected  for  back  radia¬ 
tion  from  the  model  surface. 

The  heats  of  ablation  for  teflon  and  phenolic  nylon  in  each  test  gas  are 
shown  in  Figures  5.8  and  5.  9  as  a  frnction  of  (h^  -  h^).  The  range  of  data 
is  shown  by  a  bar  in  some  cases  where  variation  in  values  were  obtained  for 
the  different  times  of  test  at  which  the  data  was  evaluated. 

The  heats  of  ablation  for  teflon  in  air  show  reasonable  agreement  with 
the  theory  of  Scala  (Ref.  25),  Good  agreement  was  obtained  at  low  enthalpy 
levels  by  Diaconis,  et  al  (Ref.  26  ). 

The  heat  of  ablation  values  for  phenolic  nylon  exhibit  a  reasonable 
trend  with  respect  to  the  approximate  theoretical  curve  shown.  It  is  pointed 
out  that  the  heat  of  ablation  for  a  charring  material  depends  upon  the  distri¬ 
bution  of  the  various  parts  of  the  total  heat  input.  Thus,  the  heat  of  ablation 
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can  vary  depending  upon  the  char  thickness  and  its  variation  with  time,  the 
virgin  material  boundary  wall  temperature,  the  rate  of  oxidation,  the  mole¬ 
cular  weight  of  the  transpired  gases  and  other  factors  that  are  difficult  to 
include  accurately  in  a  theoretical  description.  The  theory  shown  is  for  an 
"idealized"  charring  ablator  and  is  not  meant  to  represent  an  accurate 
theory  for  the  teat  case. 

5.  5  Conclusions 

The  stagnation  point  steady  state  ablation  performance  of  two  general 
types  of  ablators  were  evaluated  in  a  high  purity-high  enthalpy  free  jet  arc 
facility  in  air  and  nitrogen-carbon  dioxide  gas  mixtures.  The  mixtures  are 
assumed  to  represent  possible  Mav  and  Venus  atmospheres.  The  ratio  of 
radiative  to  convective  heat  transfer  was  small  and  therefore  the  results 
do  not  reflect  material  performance  in  a  highly  radiating  environment. 

The  heat  of  ablation  results  suggest  that  ablation  is  more  severe  in 
air  and  the  25%  CO2  mixture  than  in  N2  and  the  9%  CO^  mixture;  however, 
the  maximum  spread  between  data  points  is  only  about  10%  for  phenolic 
nylon  and  20%  for  teflon.  Thus,  it  is  concluded  that  the  steady  state  abla¬ 
tion  performance  does  not  differ  grossly  between  air  and  the  simulated 
planetary  atmospheres  at  least  for  the  materials  and  test  conditions  investi¬ 
gated. 

The  test  results  for  graphite  were  found  to  be  in  error  because  of  the 
large  mass  losses  experienced  during  cooldown  after  the  test  flow  was  ter¬ 
minated.  However,  the  mass  losses  were  still  less  than  those  found  for 
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phenolic  nylon,  Alto,  the  surface  temperatures  for  graphite  measured 
during  the  test  flow  were  significantly  higher  than  those  for  phenolic  nylon. 
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TABLE  5.  1 

ARC  TEST  CONDITIONS 


Air 

Nz 

10%  CO, 

25%  CO; 

High  Enthalpy  Test  Conditions 

Pg  (psig) 

3.0 

2.8 

2.9 

3.  1 

h/RT^ 

364 

415 

400 

384 

ho/RTo 

3.9 

3.9 

-11.  5 

-35.  5 

h3/RT^ 

368 

419 

388 

348 

h^/RTo 

(T  = 

'  w 

1500°R) 

11.  4 

11.4 

-4.2 

-28.  3 

I's  ■  Nv 

(Tw  = 

1500OR) 

357 

408 

393 

377 

RT^ 

o 

"ot 

(T  = 
'  w 

1500°R) 

2480 

2900 

2900 

2920 

(T  = 

5400°R) 

52 

45.6 

40.  3 

26.8 

hg  -  h^ 

(Tw  = 

5400°R) 

316 

37  3 

348 

322 

RTo 

"ot 

(T  = 

'  w 

5400°R) 

2260 

2610 

2540 

2480 

Low  Enthalpy  Test  Conditions 

Pg  (P«g) 

2.  3 

2.4 

2.  5 

h/RT 

o 

286 

306 

306 

s  o 

290 

295 

270 

I's  ■  Nv 

(T 

\ 

=  1500°R) 

279 

299 

299 

RTo 

qot 

(T 
'  w 

=  1500°R) 

1940 

2120 

2240 

I's  "  IV 

(Tw^ 

=  5400°R) 

244 

254 

244 

RTo 

qot 

(Tw 

=  5400°R) 

1740 

1840 

1840 

Ro  =  1/4" 

RT  =  33.86  BTU/LB 
o 

■  (T  ,  =  5400°R)  =  440  BTU/LB  (Back  Radiation  from  surface) 

Mrad  ^ 
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TABLE  5.  2 
TEST  RESULTS 


Material 

Gas 

h./RT^ 

Time 

(sec) 

(gms) 

t  char 
(in.) 

Phenolic 

Air 

368 

2.4 

.  0404 

.  029 

•  • 

Nylon 

3 

.  0505 

.032 

.67 

4 

.  0522 

.  027 

.  56 

5 

.0615 

.034 

.  47 

6.  1 

.  0745 

.  040 

.25 

419 

2 

.  0351 

.050 

1 

3 

.  0382 

.. 

6 

.  0730 

10%  CO 2 

388 

2 

.  0303 

.052 

.  80 

4 

.  0463 

.055 

.  64 

5 

.  0628 

.047 

.  56 

6 

.  0652 

.  066 

.  47 

25%  CO2 

348 

2.5 

.0434 

.  071 

.76 

4 

.  0613 

.053 

.  58 

6 

.0765 

.  026 

.42 

Teflon 

Air 

368 

2.3 

.  1358 

.94 

3.9 

.  2005 

.89 

6.2 

.  2726 

.73 

419 

2 

.  0865 

.97 

La 

3 

.  1218 

.94 

6 

.  2361 

.86 

10%  CO2 

388 

2 

.  0836 

.97 

4 

.2063 

.  91 

6 

.  3166 

.84 

25%  CO, 

348 

2.2 

.0884 

.91 

C 

4,  2 

.  1872 

.80 

5.8 

.2955 

.71 

Phenolic 

^2 

290 

2.  2 

.0276 

.  072 

.86 

Nylon 

4 

.  0458 

.050 

.67 

10%  CO2 

295 

2 

.  0293 

.  055 

.  80 

4 

.0464 

.  073 

.76 

25%  CO2 

270 

2 

.  0284 

.044 

.84 

4.  1 

.  0494 

.052 

.80 

Teflon 

Nz 

290 

2 

.  1858 

.97 

4 

.  187 

.97 

10%  CO  2 

295 

2.  1 

.  0834 

.97 

4.3 

.  1728 

.97 

25%  CO2 

270 

2 

.  0873 

1.  00 

4.  1 

.  1887 

1.00 
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THROTTLING  VALVE 


Figure  5.  2  Model  Testing  -  Tandem  Gerdien  Free  Jet  Arc 


Fij^urc  5.  5  Arc  H»Mter  -  Tanclem  Gcrdicn  Free  Jet 
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Figure  5.  4  Total  Enthalpy  Calorimeter 
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igure  5.  5  Typical  Model 
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Figure  5.  6  Typical  Model  Specimens 


Figure  5.  7a  Mass  Loss  of  Ablators 
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Figure  5.  7c  Mass  Lota  of  Ablators 
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Figure  5.  8  Heats  of  Ablation  for  Teflon 
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Figure  5.  9  Heats  of  Ablation  for  Phenolic  Nylon 


II1-6  VENUS  IONOSPHERIC  ATTENUATION* 


6.  1  Introduction 

The  object  of  this  part  of  the  study  program  was  to  determine  the 
expected  amount  of  attenuation  of  a  radio -frequency  signal  passing  through 
the  ionosphere  of  Venus.  The  attenuation  of  a  signal  in  the  ionosphere  de¬ 
pends  on  the  electron  concentration  and  the  collision  frequency,  which  in 
turn  depend  quite  critically  on  the  atmospheric  composition.  The  composi¬ 
tion  data  used  for  the  three  models  are  those  given  in  Section  I  of  this  report. 
The  attenuation  has  been  calculated  for  transmission  frequencies  of  100  me 
and  500  me. 

Initial  considerations  of  the  atmospheric  models  provided  by  JPL  in¬ 
dicated  that  losses  in  the  ionosphere  would  be  sufficiently  small  that  mathe¬ 
matically  convenient  assumptions  might  yield  results  of  acceptable  accuracy 
for  the  present  study.  The  results  appear  to  justify  these  assumptions.  The 
two  basic  assumptions  are  that  the  propagation  can  be  treated  in  a  piece - 
wise-linear  manner  and  that  reflection  at  the  boundaries  can  be  neglected. 

The  resulting  expression  is; 


P  =  P  ?  “"‘i’  '‘'“’i 
o 

where  P^  is  the  power  incident  on  the  ionosphere  andQ(z) 
coefficient,  given  by  (Ref.  28); 


(6-1) 

is  the  attenuation 


a(z)  =  (fp/f)^i/(z)  (1 


2  2-1 


-1 

c 


(6-2) 


where. 


f  =  e  V  N  /irm 
P  e 

■fPrincipal  Investigators:  S.  Visvanathan  and  W.  C.  King 


(6-3) 
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is  the  critical  frequency 

is  the  electron  concentration,  given  by  JPL  as  10^/ cc. 
f is  the  signal  frequency 
c  is  the  velocity  of  light 

V{z)  is  the  collision  frequency  of  the  electrons 

6.2  Collision  Frequencies  for  the  Three  Atmospheric  Models 
The  collision  frequency  is  given  by  (Ref.  29); 

u{z)  =  N(z)a^^^V  (6.4) 

where, 

N  (z)  is  the  density  of  the  electrically  neutral  particles  at  altitude  z 
is  the  effective  cross-section  of  the  electrons 
T  is  the  temperature  of  the  electrons 
V  is  the  mean  velocity  of  the  electrons,  given  by 


The  effective  collision  cross  section  for  the  Venus  model  atmospheres 
is  subject  to  considerable  uncertainties.  As  observed  from  experimental 
work  reported  in  Massey  &  Burhop  (Ref.  29),  and  also  in  Brown  (Ref.  30), 
experimental  data  in  the  energy  ranges  of  interest  (considerably  less  than 
one  electron  volt)  is  sufficiently  scarce  as  to  be  unreliable.  Extrapolation 
of  experimental  data  permits  a  reasonable  choice  within  a  range  of  about 
an  order  of  magnitude. 

For  the  purposes  of  this  study,  the  effective  cross  sections  for  the 
three  model  atmospheres  are  taken  as  follows; 
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7.  64  X  10“^^  cm^ 

7.  2  X  10”  cm^ 

8.73  X  10*^^  cm^ 

The  collision  frequencies  were  computed  from  Equation  (6-4),  using 
these  values  of  the  cross  section  and  the  corresponding  values  for  the  neu¬ 
tral  density  N  and  the  thermal  velocity.  The  results  of  these  calculations 
are  given  in  Table  6.  1. 


Best  Model; 

•^eff  = 

Ebctreme  I; 

Q 

11 

Extreme  II; 

'^eff  = 

121 


Height  in  km 
110 
122 
133 
144 
178 
235 
292 

Height  in  km 

no 

122 

133 

144 

178 

235 


Height  in  km 
100 
150 
200 
250 
300 


Table  6.  1 
Best  Atmosphere 

Collision  Frequency  no/ sec 

7.  52  X  105 
1.  31  X  10® 

3.62  X  10^ 

1.  19  X  10^ 

8,  12  X  10^ 

8.  16  X  10^ 

0.  11 


Extreme  I 


Extreme  II 


122 


Collision  Frequency  no/ sec 
9.  58  X  10® 

2.  94  X  10® 

3.  67  X  10^ 

6.44  X  10® 

7.  70  X  10^ 

2.  73  X  10^ 


Collision  Frequency  no/ sec 

1.69  X  10® 

0.  6 

2.  28  X  lO-** 

6xl0-'0 
3  X  lO"^** 


6. 3  Attenuation  in  the  Venusian  Ionosphere 


The  critical  frequency  for  the  Venusian  ionosphere,  taking  the  electron 


concentration  as 

10^/ cc  is 

f  = 

P 

8.  984  me 

(6-6) 

The  attenuation  0((z)  has  been  computed  for 

the  three  models  using  equation 

(6-2)  and  the  above  data.  The  results  are 

presented  in  Table  6.  2  below  for 

the  signal  frequencies  f  =  100  me  and  500  me. 

Table  6.  2 

Best  Atmosphere 

Height  in  km 

Attenuation  in  cm"^ 

Attenuation  in  em"^ 

(f  =  100  me) 

(f  =  500  me) 

no 

1.023  X  10*'^ 

3.76  X  10-9 

122 

1.78  X  10-® 

6.  56  X  10“ 

133 

4.  92  X  10‘9 

1.81  X  10"^° 

144 

1.62  X  10"^ 

5.  94  X  10-^^ 

178 

1. 10  X  10*^0 

4.06  X  10“^^ 

235 

1 

0 

X 

4. 08  X  10"^^ 

292 

1.  50  X  lO'^*^ 

Extreme  I 

5.  50  X  10-^^ 

Height  in  km 

Attenuation  in  em'^ 

Attenuation  in  cm"^ 

(f  =  100  me) 

(f  =  50  me) 

no 

1. 30  X  10"^ 

4.79  X  10-9 

122 

4.00  X  10"® 

1.47  X  10-9 

133 

1.40  X  10‘® 

5. 14  X  10-^9 

144 

4.  50  X  lO'*^ 

1.84  X  10-10 

178 

8.76  X  lO'^O 

3. 22  X  10-11 

235 

1.05  X  10"^0 

3.85  X  10-1^ 

292 

3.71  X  10"^^ 

1.  37  X  10-1^ 
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Extreme  II 


Height  in  km 

Attenuation  in  cm  '  ^ 

Attenuation  in  cm 

(f  =  100  me) 

(f  =  500  me) 

100 

2.30  X  10’® 

8,45  X  10-^0 

150 

8. 16  X  10"^^ 

3,  08  X  10'^5 

200 

3.82  X  lO-^"^ 

1. 15  X  10-^® 

The  total  attenuation  was  obtained  by  summing  the  contributions  from 
the  altitude 'increments: 

''total  =  \ 

The  calculated  attenuations  in  decibels  are  given  in  Table  6.  3  for  the 
three  atmospheric  models. 


Table  6.  3 

Total 

Total 

Model 

Path  Length 

Attenuation  at 

Attenuation  at 

100  me 

500  me 

Best 

182  km 

0.80  db 

0.  02  db 

Extreme  I 

182  km 

0.  52  db 

0.017  db 

Extreme  II 

100  km 

1.  35  db 

0.  05  db 

These  attenuations,  of  course,  are  for  a  vertical  propagation  path  through 
the  ionosphere.  The  attenuation  vs.  path  length  and  altitude  is  shown  in 
Figures  6.  1  and  6.  2. 

Table  6.  3  shows  that  there  is  no  problem  in  communicating  through 
the  ionosphere  of  Venus.  The  Extreme  II  model  yields  the  pessimistic  esti¬ 
mate,  and  Extreme  I  yields  the  optimistic  estimate.  It  should  be  pointed 
out,  however,  that  it  is  not  known  exactly  where  the  ionosphere  of  Venus 
begins;  it  has  been  estimated  that  it  starts  around  100  km.  For  a  true  eval¬ 
uation  of  the  attenuation,  one  should  know  the  variation  of  the  electron  con- 
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cent  ration  with  altitude,  particularly  around  100  km,  where  the  collision 
frequency  gets  higher;  this  aspect  needs  to  be  investigated  more  carefully 
when  better  knowledge  of  the  composition  of  the  Venusian  atmosphere  be¬ 


comes  available. 
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Figure  6.  1.  Venus  Ionospheric  Attenuation  vs.  Path  Length  at  f  =  100  me. 
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Figure  6.  2.  Venus  Ionospheric  Attenuation  vs.  Path  Length  at  f  =  500  me. 


IV.  RECOMMENDED  FOLLOW-ON  RESEARCH  AND  DEVELOPMENT 


PROGRAM 

Though  not  all  of  the  planetary  entry  problems  could  be  solved  in  a 
single  short  program,  this  program  has  accomplished  the  following:  First, 
it  has  achieved  the  majority  of  the  objectives  established  for  the  program 
(see  section  II).  Second,  it  has  indicated  additional  real  and  potential  prob¬ 
lem  areas.  Third,  it  has  established  techniques  and  procedures  for  future 
investigation  of  these  problems.  The  following  recommended  programs  do 
not  include  all  of  the  recognized  problem  areas,  but  rather  represent  those 
that  should  be  pursued  as  the  recommended  follow-on  studies  of  this  program. 

IV.  1  EQUILIBRIUM  THERMOCHEMISTRY 

The  Rand  method  program  should  undergo  extensive  check-out  and 
verification  with  the  results  compared  to  those  of  the  Brinkly  method  pro¬ 
gram.  Subsequent  to  the  verification,  equilibrium  thermochemistry  tables 
should  be  prepared  for  the  mixtures  where  greater  accuracy  is  desired 
than  is  available  from  the  Brinkly  method.  Upon  establishment  of  the  Mars 
and  Venus  atmospheres,  tables  and  Mollier  type  charts  should  be  prepared 
for  these  atmospheres.  Charts  of  this  nature  will  be  utilized  as  major  ref¬ 
erences  similar  to  the  AVCO  chart  for  air  which  is  prevalently  used  by  the 
industry  today. 

IV-2  DISSOCIATION  RATE  STUDIES 

With  the  completion  of  the  study  of  the  dissociation  of  nitrogen  behind 
combustion  generated  shock  waves,  an  investigation  should  be  continued 
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to  determine  the  influence  on  these  results  that  would  be  caused  by  the  pre¬ 
sence  of  small  quantities  of  CO2.  This  would  provide  information  pertinent 
to  the  atmosphere  of  Mars.  Additional  investigations  should  be  undertaken 
of  the  CO  system,  then  of  the  CO2  system  and  then,  at  a  later  date,  an  exten¬ 
sive  investigation  should  be  made  of  the  combined  carbon  dioxide -nitrogen 
system. 

At  the  contemplated  entry  velocities  of  40,  000  feet  per  second  into 
the  atmosphere  of  Venus  there  should  be  enough  thermal  energy  available  to 
partially  dissociate  CO,  which  is  one  of  the  decomposition  products  of  CO2. 
Since  we  have  recently  completed  a  shock  tube  study  of  the  dissociation  of 
nitrogen  (Ref.  6)  (which  is  the  chief  component  of  the  atmospheres  of  Mars 
and  Venus),  it  is  proposed  to  study  the  thermal  dissociation  of  CO  behind 
strong  shocks. 

While  it  is  generally  believed  that  the  ambient  Martian  atmosphere 
contains  chiefly  nitrogen  and  some  carbon  dioxide,  a  recent  paper  by  Opik 
(Ref.  31)  has  suggested  that  the  CO^  on  Mars  may  be  decomposed  into  CO 
and  O  by  ultraviolet  absorption.  However,  the  work  of  Mahon  (Ref.  32)  tends 
to  contradict  this  suggestion,  indicating  that  photolysis  of  CO2  leads  to  the 
formation  of  ozone,  which  would  prevent  further  break  up  of  the  CO2  by  solar 
radiation.  In  any  event, CO  will  be  present  in  the  flow  surrounding  a  vehicle 
and  will  probably  be  partially  dissociated. 

The  2. inch  square  shock  tube  which  utilizes  combustion  driving  de¬ 
scribed  in  section  III-2  is  well  suited  to  the  task  of  studying  the  thermal  dis- 


129 


Bociation  of  CO.  In  the  case  of  the  disBociation  of  the  COithe  carbon  atom' s 


first  ionization  potential  lies  only  slightly  above  the  dissociation  energy  of 
the  CO.  Therefore,  there  is  the  interesting  possibility  that  some  of  the  car¬ 
bon  atoms  may  be  ionized.  If  there  is  enough  energy  available  to  form  ct  then 
there  will  be  enough  to  excite  the  O  atoms.  The  oxygen  atom  has  a  line  that 
is  particularly  well  suited  for  Stark  broadening  measurements.  A  method 
has  been  formulated  which  permits  the  calculation  of  the  Stark  broadening. 

With  combustion  driving  it  will  be  possible  to  dissociate  up  to  10%  of  the  CO 
behind  shocks  in  pure  CO  and  up  to  70%  in  CO-Ar  mixtures.  On  the  basis  of 
the  recent  nitrogen  (Ref.  6)  study  it  is  likely  that  the  temperature  may  fall 
by  2,  000°K  in  the  reaction  zone.  Such  a  temperature  drop  will  probably  pro¬ 
hibit  much  C**  formation.  If  very  little  C**  is  formed, then  the  dissociation 
process  is  most  likely  bimolecular.  The  program  written  for  the  IBM  1620 
in  the  nitrogen  experiment  can  be  readily  applied  to  CO. 

An  additional  study  that  will  produce  early  results  with  a  minimum  of 
expended  effort  is  an  investigation  of  oxygen.  The  nitrogen  results  are  in 
disagreement  with  previous  results  for  nitrogen  (Ref.  7)  and  differ  from  the 
results  on  oxygen  (Refs.  33  and  34).  It  is  suggested  from  this  that  the  pre¬ 
vious  results  (Refs.  33  and  34)  inferring  that  O  is  some  3-5  times  more  effi¬ 
cient  than  O2  in  dissociating  O2  may  be  in  error. 

Knowing  the  relative  efficiency  of  O  in  dissociating  O2  is  of  applied 
interest  for  space  vehicles.  It  has  been  shown  (Ref.  35)  that,  of  all  the  rates 
in  the  air  system,  electron  production  is  most  sensitive  to  the  rate  of  disso- 
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elation  of  oxygen.  In  fact,  a  change  by  a  factor  of  10  in  the  oxygen  dissocia¬ 
tion  rate  leads  to  a  change  in  the  rate  of  production  of  electrons  by  200.  At 
stagnation  temperatures  of  6,  400°K  the  O2  almost  completely  dissociates; 
hence,  the  oxygen  molecules  are  dissociated  principally  by  collisions  with 
oxygen  atoms  and  nitrogen  molecules.  It  is  clear  that  changing  the  rate  of 
dissociation  of  O2  through  collisions  with  O  by  a  factor  of  3  may  lead  to  a 
change  in  the  over -all  rate  of  production  of  electrons.  In  vehicles  flying  at 
high  altitudes  a  change  by  an  order  of  magnitude  in  the  rate  of  production  of 
electrons  could  move  the  point  of  equilibrium  electron  concentration  a  size¬ 
able  distance  relative  to  the  vehicle.  Such  a  result  could  have  a  profound 
influence  on  the  choice  of  antenna  location,  radio  frequency  and  power  re¬ 
quirements. 

It  is  not  proposed  here  to  do  a  complete  study  of  oxygen  dissociation, 
but  rather  to  make  a  careful  check  to  see  if  the  reported  rates  (Refs.  33  and 
34)  will  properly  fit  the  experimental  density  profiles  under  conditions  of 
good  experimental  accuracy.  Such  data  should  clarify  the  question  of  the 
reported  anomalously  high  efficiency  of  the  O  atom.  If  the  previously  re¬ 
ported  values  are  not  found  to  be  sufficiently  accurate  to  predict  the  density 
profile,  then  a  more  exhaustive  study  would  be  in  order. 

IV -3  BLUNT -BODY  STAGNATION-POINT  HEAT  TRANSFER  - 
EXPERIMENTAL  STUDY 

Two  problems  have  not  been  investigated  fully  in  the  present  study 
which  will  have  bearing  upon  the  stagnation  region  heating  problem  of  a  blunt 
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body  entering  the  atmosphere  of  Venus  (or  Mars).  First,  the  radiative  heat 
transfer  may  be  a  large  contributor  to  the  total  heat  transfer.  Second,  the 
flow  behind  the  bow  shock  wave  and  in  the  shock  layer  may  not  be  in  equili¬ 
brium  for  important  ranges  of  flight  conditions.  This  could  add  significantly 
to  the  radiative  heat  transfer,  if  there  are  radiation  overshoots  in  the  non¬ 
equilibrium  regions,  and  may  appreciably  affect  the  convective  or  boundary 
layer  heat  transfer  as  well.  The  influence  of  small  concentrations  of  water 
vapor  on  the  non-equilibrium  characteristics  of  simulated  Venus  atmosphere 
should  be  considered  as  an  important  part  of  this  program. 

Therefore,  the  following  two  study  are  recommended.  First,  an 
experimental  program  should  be  undertaken  to  measure  the  equilibrium  and 
non-equilibrium  heat  transfer  rates  from  shocked  gases  over  a  range  of  mix¬ 
tures  (including  small  concentrations  of  water  vapor  or  hydrogen  to  simulate 
water)  that  simulates  the  atmospheres  of  Venus  and  Mars.  The  extent  of  non¬ 
equilibrium  regions,  the  properties  and  concentrations  of  gases  within  them, 
and  their  dependence  upon  flight  conditions  would  be  determined  as  part  of 
this  program.  Emphasis  would  be  placed  upon  total  radiation  measurements 
(with  a  cavity  gage),  but  time  resolved  spectroscopic  observations  would  also 
be  made  to  identify  the  radiating  species  and  their  importance  in  both  the 
equilibrium  and  non-equilibrium  regimes.  These  results  would  be  compared 
with  existing  theoretical  predictions.  The  effects  of  body  flow  field  property 
gradients  would  also  be  studied.  Second,  the  present  convective  heat  trans¬ 
fer  measurements  should  be  extended  to  wider  ranges  of  variables.  Also, 
the  effects  of  boundary  layer  and  shock  layer  non-equilibrium,  as  identified 
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by  the  radiation  study,  on  the  convective  heating  should  be  determined. 

The  two  recommended  studies  should  be  conducted  in  a  high  perform¬ 
ance  shock  tube  that  is  capable  of  producing  test  flows  in  simulated  atmo¬ 
spheric  mixtures  over  a  range  of  pressure  levels  and  to  velocities  beyond 
40,  000  ft/ sec.  The  new  GE-MSD  6"  hypervelocity  shock  tube  which  is  driven 
by  arc -heated  helium  (300  KJ  power  supply)  is  such  a  facility. 

IV-4  BLUNT-BODY  LAMINAR  AND  TURBULENT  HEAT  TRANSFER- 
EXPERIMENTAL  STUDY 

Ihe  radiative  and  convective  heat  transfer  to  surfaces  away  from  the 
stagnation  region  of  a  body  shotild  be  investigated  for  conditions  and  gas  mix¬ 
tures  simulating  Venus  (and  Mars)  entry  flight.  Convective  heat  transfer 
distributions  should  be  measured  as  a  function  of  body  shape  and  size  and 
for  a  range  of  simulated  altitudes  and  velocities  to  greater  than  40,  000  ft/ 
sec.  Model  pressure  levels  should  be  varied  at  an  appropriately  high  enthalpy- 
level  to  cover  the  range  fromfxlly  laminar  to  turbulent  body  boundary  layer. 

The  date  should  be  compared  with  existing  theoretical  predictions  and  plotted 
against  appropriate  correlation  parameters.  Thus,  the  suitability  of  current¬ 
ly  usefvil  analytical  and  empirical  techniques  for  the  gases  in  question  at  the 
high  stagnation  enthalpy  levels  of  interest  would  be  assessed.  A  possibly 
important  result  of  this  work  would  be  the  production  of  new  boundary  layer 
transition  data  v^ich  would  be  compared  with  results  now  available  for  lower 
flight  velocity  simulation.  Because  of  the  high  pressures  and  large  model 
sizes  available  with  the  MSD-SSL  6"  shock  tube,  it  should  be  possible  to 
produce  a  turbulent  boundary  layer  on  the  forward  portions  of  a  blunt  shape 
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(that  if,  in  the  regions  of  high  free  stream  temperature  and  ionization  levels) 
so  that  the  effects  of  pressure  gradient  or  turbulent  heating  and  transition 
can  be  studied. 

Total  radiative  heating  to  surfaces  away  from  the  stagnation  region  of 
a  body  should  be  measured  to  determine  the  extent  of  the  radiative  heating 
problem.  This  would  be  done  on  bodies  of  various  shapes  and  sizes  using 
optically  measured  shock  layer  shapes  and  surface  total  heat  transfer  gages 
(cavity  gage).  These  experiments  would  be  conducted  in  high  enthalpy  flows 
in  which  the  shock  layer  is  dominated  by  equilibrium  radiation.  Approximate 
analytical  methods  would  be  developed  using  previously  determined  non-equi¬ 
librium  radiative  data  for  normal  shock  waves  to  construct  the  complete 
radiative  heating  distributions  for  bodies  in  atmospheric  flight  over  the 
desired  ranges  of  altitudes  and  velocities. 

The  recommended  experimental  study  would  be  conducted  in  the  GE- 
MSD  6"  hypervelocity  shock  tube.  This  facility  is  driven  by  arc-heated 
helium  (300  KW  power  supply).  This  tube,  in  addition  to  being  capable  of 
producing  flight  velocity  simulation  to  greater  than  40,  000  ft/ sec,  can  also 
produce  and  contain  high  stagnation  pressure  levels  which  are  required  for 
turbulent  boundary  layer  studies.  (The  order  of  100  atmospheres  at  greater 
than  30,  000  ft/ sec  simulated  flight  velocity  is  possible  with  the  present  tube 
design  and  power  supply. ) 
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IV-5  AEROTHERMODYNAMIC  PERFORMANCE  OF  HEAT  PRO¬ 


TECTION  MATERIALS  -  EXPERIMENTAL  STUDY 

5.  1  Introduction 

The  heat  protection  materials  study  conducted  for  the  present  program 
considered  the  relative  performance  of  several  materials  representative  of 
various  classes  of  ablators  exposed  to  high  levels  of  convective  heating  in 
the  subject  simulated  environments.  Trajectory  analyses  performed  in  the 
recent  Venus-Mars  Capsule  Study,  conducted  by  the  Missile  and  Space  Divi¬ 
sion  of  the  General  Electric  Co.  under  contract  to  the  NASA  Jet  Propulsion 
Lab,  clearly  point  out,  as  expected,  that  the  extremely  high  velocities  asso¬ 
ciated  with  the  Venus-Mars  re-entry  problem  will  impose  an  appreciable 
radiative  heating  load  upon  the  stagnation  region  of  the  vehicle;  this  can  be 
a  dominating  factor  in  the  heat  protection  design  (Ref.  18).  This  facet  of  the 
materials  evaluation  was  not  examined  in  the  above  study.  In  addition,  con¬ 
sideration  of  the  skirt  and  aft  portions  of  such  a  vehicle,  where  heat  trans¬ 
fer  is  an  order  of  magnitude  or  more  lower  than  stagnation  point  values, 
can  also  impose  thermal  protection  problems  unlike  those  examined  pre¬ 
viously.  Material  performance  in  this  latter  case  can  be  classified  as  non¬ 
steady  ablation  in  contrast  to  the  steady  state  ablation  phenomenon  exempli¬ 
fied  by  exposure  to  high  heat  transfer  environments  (Ref.  36).  Whereas  for 
steady  state  ablation  the  damage  of  heat  penetration  of  the  material  proceeds 
at  the  same  rate  as  the  surface  recedes,  for  so-called  non-steady  ablation 
the  heat  penetrates  into  the  material  at  a  faster  rate  than  the  surface  recedes. 
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Phyaically  for  ateady  atate  ablation  there  reanlta  a  damage  layer  that  ia  thin 
but  movea  very  rapidly  into  the  material,  a  aituation  in  which  primarily 
maaa  losa  rather  than  inaulation  capability  ia  the  critical  ablation  property. 

On  the  contrary,  in  non-ateady  ablation  the  damage  layer  ia  rather  thick 
(aee  Figure  IV-5.  1)  and  aurface  recesaion  ia  quite  alow, making  the  critical 
criterion  for  evaluation  of  a  material  the  amount  of  thicknesa  or  depth  re¬ 
quired  to  maintain  a  certain  level  of  backwall  ahield  temperature.  Both 
aapecta  of  the  heat  protection  problem,  this  non-steady  ablation  phenomenon 
and  the  serious  radiation  heat  transfer  problem  previously  mentioned,  must 
be  thoroughly  examined  before  an  optimum  heat  protection  design  is  attempted 
for  a  Venus  or  Mars  entry  capsule.  Such  a  study  is  described  in  the  ensuing 
discussion. 

Before  embarking  on  any  proposed  materials  evaluation  program  it 
is  advisable  to  define  the  free  flight  environment  of  interest.  Shown  in  Figures 
IV-5.  2  and  IV>5.  5  are  typical  stagnation  point  heat  transfer  profiles  deter¬ 
mined  for  the  various  entry  trajectories  considered  in  the  Venus-Mars  cap¬ 
sule  study  performed  by  MSD  in  the  recent  study  for  JPL  (Ref.  18).  Included 
are  calculations  for  various  model  atmospheres  and  vehicle  shapes.  A 
perusal  of  these  data  show  that  over  the  range  of  entry  parameters  considered 
and  indicated  on  the  various  figures,  peak  radiative  heating  rates  might  vary 
from  approximately  100  BTU/ft  -sec  to  30,000  BTU/ft  -sec  for  a  sphere- 
cone  body  1.  5  ft,  in  radius, while  the  convective  values  are  bracketed  between 
2,000  and  20,000  BTU/ft^-sec.  Further,  the  total  heating  environment 
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(combined  radiative  and  convective)  liee  between  approximately  30,  000  and 
70,  000  BTU/ft  for  a  1  1/2  ft.  rad.  body.  The  magnitude  of  these  heating 
rates  is  extremely  large,  and  the  design  problem  a  serious  one  yet  not  much 
work  has  been  performed  in  ground-based  facilities  to  evaluate  materials 
under  such  conditions.  The  lack  of  extremely  high  enthalpy/ hs  700 

Wo 

material  test  facilities  has  precluded  the  attainment  of  such  values,  although 
short  time  aerodynamic  shock  tunnel  facilities  with  such  capabilities  are 
presently  in  operation  at  MSD  and  elsewhere  (Refs.  37  and  38).  As  is  well 
known,  long  time  electric  arc  facilities  are  required  to  provide  the  necessary 
simulation  environments  for  the  evaluation  of  ablation  materials,  at  least 
for  the  establishment  of  severe  convective  heating  test  conditions.  Simula¬ 
tion  of  large  radiative  heating  loads  in  the  laboratory  is  not  as  far  advanced, 
due  primarily  to  only  the  recent  requirement  for  such  environments  for 
study  purposes.  In  addition,  since  the  radiative  loads  attainable  are  directly 
proportional  to  the  volume  of  the  radiating  species  (say  test  gas),  the  use  of 
small  test  models  in  laboratory  test  facilities  has  further  complicated  the 
problem  of  the  attainment  of  large  radiative  heat  transfer  values  on  test 
specimens. 

Recently  at  the  SSL  a  test  technique  has  been  developed  whereby  it  is 
feasible  to  conduct  studies  in  controlled  radiative  and  convective  heating  en¬ 
vironments  using  the  tandem -Gerdien  arc  heater  developed  at  this  laboratory 
(Ref.  39).  Unlike  most  approaches  to  this  problem  which  are  handled  by  a 
superposition  of  a  separate  radiative  source  upon  a  plasma  jet  flow,  this 
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technique  makes  use  of  a  single  electric  arc  source  for  both  the  radiative 

and  convective  heating  loads.  A  salient  advantage  of  this  arrangement  is  that 

the  source  of  the  radiation  is  high  temperature  air  resulting  in  a  spectral 

energy  distribution  similar  to  that  which  would  be  experienced  by  entering 

vehicles  at  high  velocity.  The  technique  has  been  evaluated  and  correlated 

2 

at  low  radiant  intensities  (^200  Btu/ ft  sec)  under  conditions  whereby  the 
gas  enthalpy  (at  one  atmosphere  pressure)  measured  at  the  model  specimen 
location  was  found  to  be  around  450  RT^.  Subsequent  studies  have  resulted 
in  gas  stagnation  enthalpy  measurements  in  excess  of  800  RT^  with  certain 
test  configuration  modifications.  Currently  calibration  measurements  are 
being  made  at  environmental  pressure  levels  of  five  atmospheres  which, 
when  coupled  with  the  higher  enthalpy  test  flows,  should  provide  a  test  radia- 
tion  intensity  in  excess  of  1000  Btu/ft  sec  at  the  specimen  surface.  Simultan¬ 
eously  with  this  large  radiant  source  a  controllable  (depending  upon  amount 
of  by-passed  flow)  convective  heating  load  from  the  flowing  high  enthalpy 
gases  can  also  be  superimposed  on  the  specimen.  In  the  following  section 
a  more  detailed  description  of  the  test  technique  is  presented. 

5.  2  The  Tandem-Gerdien  Radiation  Arc 

The  product  of  considerable  development  effort,  the  Space  Sciences 
Laboratory  tandem-Gerdien  arc  heater  produces  a  high  enthalpy  air  flow 
with  an  insignificant  trace  (^  100  ppm)  of  carbon  contaminant  from  the  elec¬ 
trodes  for  operating  times  in  excess  of  1  5  minutes.  This  configuration  has 
been  operated  with  various  test  facilities  (wind  tunnel  and  free  jets)  in  the 
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Laboratory  consistently  at  air  stagnation  enthalpy  values  (^s/RT^)  in  excess 
of  400;  however,  to  produce  high  radiant  intensity  levels  for  material  testing, 
calculations  show  that  enthalpies  in  excess  of  700  RT^  are  desired  (Ref.  39). 
Data  obtained  during  the  developmental  period  of  the  arc  heater  indicated 
that  the  necessary  enthalpy  requirement  could  be  achieved  if  the  test  gas 
specimen  were  extracted  from  the  arc  column  where  gas  temperatures  in 
excess  of  15,  000°K  had  been  measured  by  McGinn  (Ref.  40).  After  consider¬ 
able  study  and  design  evaluation  the  configuration  shown  in  Figure  IV-5.  6  was 
constructed.  A  photograph  of  the  facility  is  shown  in  Figure  IV-5.  7.  The 
major  alteration  involved  the  plenum  chamber  whose  basic  shape  was  changed 
from  a  cylindrical  section  to  a  tee  section;  the  diameter  of  all  three  legs  of 
the  tee  is  equal  to  that  of  the  restrictor  plate  orifices.  With  this  arrange¬ 
ment  the  hot  gas  that  normally  mixed  in  the  larger  plenum  in  the  presence 
of  the  arc  column  is  further  energized  by  the  very  core  of  the  column  before 
flowing  towards  the  model.  Besides  the  processed  gas  in  the  column,  calcu¬ 
lations  show  that  the  gas  flowing  towards  the  model  contains  sufficient  energy 
to  contribute  appreciably  to  the  radiant  intensity  at  the  model  surface  if  suffi¬ 
ciently  high  enthalpy  has  been  generated  by  the  arc  heater  (**s/RTq^  700). 

As  indicated  in  Figure  IV-5.  6,  the  test  configuration  is  a  shroud 
nozzle  so  that  with  a  high  enthalpy  gas  generated  in  the  plenum, large  con¬ 
vective  heat  transfer  rates  would  be  imposed  on  a  test  model  in  addition  to 
the  radiative  contributions.  (Equivalently  high  convective  values  would  also 
be  imposed  if  the  test  configuration  was  operated  as  a  free  jet). 
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To  minimize  the  convective  heating  rates  so  as  not  to  mask  the  re¬ 
sponse  of  the  materials  to  the  gas  radiation,  a  by-pass  was  incorporated  in 
the  modified  plenum  whereby  a  large  portion  of  the  hot  gas  would  be  exhausted 
overboard  just  ahead  of  the  model.  Using  this  approach  the  hot  gas  flow  re¬ 
quired  for  the  radiant  energy  transfer  is  not  disturbed, but  the  convective 
heating  can  be  controlled. 

5.  3  Proposed  Study 

To  characterize  the  performance  of  ablation  materials  in  heat  trans¬ 
fer  environments  typical  of  entry  into  the  planetary  atmospheres  of  Venus 
and  Mars,  including  the  effects  of  large  radiative  heating  loads  and  non¬ 
steady  ablation  phenomena  where  applicable,  the  following  studies  are  pro¬ 
posed.  The  work  is  conveniently  subdivided  into  two  distinct  studies  which 
are  independent  of  one  another. 

A.  Systems  Study 

An  experimental  heat  protection  systems  study  would  be  conducted, 
taking  into  account  the  broad  range  of  trajectory  parameters  that  influence 
ablation  material  behavior  and  exposing  test  specimens  in  a  controlled  manner 
to  the  necessary  test  environments  to  properly  characterize  such  materials. 
Included  in  this  study  would  be  an  investigation  of  material  degradation  in 
environments  which  provide  appreciable  radiative  heating  loads,  large  con¬ 
vective  rates  and  others  in  which  non-steady  ablation  characteristics  pre¬ 
dominate.  Although  trajectory  calculations  indicate  even  for  the  more  shallow 
entry  conditions  that  a  large  percentage  of  the  heat  transfer  over  the  rear- 
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ward  portions  of  the  vehicle  would  be  turbulent,  the  increased  heating  rates 
over  the  laminar  values  may  still  be  low  enough  to  result  in  departures  from 
steady  state  ablation.  Certainly  when  the  stagnation  point  heating  as  indi- 
cated  in  Figure  IV>5.2  is  only  100  BTU/ft  -sec,  one  could  expect  heating  rates 
on  the  skirt,  even  with  turbulent  flow. that  could  result  in  the  above-mentioned 
deviations.  It  seems  desirable,  therefore,  to  pursue  this  phenomenon. 

At  least  one  existing  material  would  be  examined,  the  emphasis  of 
the  study  being  on  the  prediction  of  the  performance  of  proven  ablation  mate¬ 
rials  (in  terrestrial  environments)  for  trajectories  typical  of  a  proposed 
mission.  In  this  respect  it  is  felt  that  phenolic  nylon  is  well  suited  to  be 
included  in  such  a  study. 

Typical  entry  trajectories  would  be  chosen  on  the  basis  of  overall 
systems  considerations;  that  is,  a  family  of  entry  paths  into  the  planetary 
atmospheres  would  be  used,  covering  a  range  of  ballistic  parameter,  entry 
angle,  entry  velocity,  etc.  For  these  results  the  respective  stagnation  point 
heating  rate  histories,  both  radiative  and  convective,  for  a  typical  vehicle 
(sphere-cone)  would  be  obtained  from  the  results  of  the  Venus-Mars  Capsule 
Study.  From  these  results  for  a  limited  number  of  cases  where  it  appears 
desirable,  heating  rate  histories  at  the  aft  section  of  the  vehicle  (skirt)  will 
also  be  computed.  In  view  of  the  uncertainties  associated  with  the  correct 
species  determination  for  the  Venus  and  Mars  atmospheres,  it  is  suggested 
to  conduct  these  studies  in  the  various  simulated  environments  currently 
being  considered  by  JPL  and  MSD.  These  are; 
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a)  3%  CO2  97%  N2 

b)  9%  CO2  91%  N2 

c)  25%  CO2  75% 

A  priori  it  would  be  decided  which  of  the  calculated  trajectories  represented 
extremes  in  expected  material  performance, thereby  providing  the  criteria  for 
simulation  in  the  test  facilities. 

The  case(s)  chosen  as  representative  of  very  severe  environments 
both  from  the  radiative  and  convective  viewpoint  would  be  simulated  in  the 
tandem -Gerdien  radiation  arc  in  which  both  types  of  heating  environments  can 
be  generated.  In  these  tests  by  varying  the  operational  test  parameters  -- 
pressure  level  (1-5  atm), gas  enthalpy  (700>^s/RT^^  400),  flow  by-pass  and 
exposure  time,  the  performance  of  the  ablation  materials  investigated  would 
be  established  for  both  radiative  and  convective  environments.  Where  possi¬ 
ble  test, exposure  would  be  extended  to  total  heating  as  inferred  by  the  chosen 
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heat  transfer  trajectories  (BTU/ft  ).  With  this  information  available,  the 
requirements  of  a  full  trajectory  simulation  test  will  be  developed  and  a 
study  performed  in  the  tandem -Gerdien  arc  wind  tunnel  (see  Figure  IV-5.  8). 
The  reference  criteria  for  these  tests  will  be  a  previously  selected  typical 
heat  transfer  history  applicable  to  the  hypothetical  vehicle  skirt  location, 
although  the  actual  test  specimen  would  be  a  blunt  axisymmetric  model  evalu¬ 
ated  at  its  stagnation  point.  The  assumptions  leading  to  the  choice  of  sim¬ 
ulated  trajectory  properties  could  thus  be  checked  and  modified  conditions 
investigated  if  necessary.  Model  parameters  would  be  material  thickness, 

mass  loss,  backwall  temperature  and  surface  temperature.  Analysis  of  the 
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results  would  include  comparison  with  calculated  material  performance  in 
both  the  theoretical  flight  trajectories  and  the  simulated  test  environments. 

The  outgrowth  of  this  study  with  appropriate,  well  characterized, 
thermal  shield  ablation  materials  should  provide  the  vehicle  designer  with 
reliable  inputs  useful  for  actual  design  of  a  thermal  shield  for  a  Venus-Mars 
entry  capsule. 

B.  Evaluation  of  Improved  Materials 

A  second  facet  of  this  proposed  study  involves  the  experimental  evalu¬ 
ation  of  new  and  improved  materials  for  use  as  thermal  shields  on  future 
Venus-Mars  entry  vehicles.  The  first  study,as  indicated,  will  consider 
currently  available  ablation  materials  which  are  well  characterized  and  have 
proven  to  perform  successfully  on  terrestial  re-entry  vehicles.  For  plane¬ 
tary  entry,  with  the  anticipated  severe  environmental  conditions,  calculations 
show  that  the  thermal  shield  can  be  an  appreciable  portion  of  the  complete 
capsule  weight  (Ref.  18).  The  use  of  better  performing,  lighter  total  weight 
ablation  materials  is,  therefore,  desirable  to  improve  mission  capabilities 
and  assure  successful  entry.  Areas  of  contemplated  material  research  leading 
to  more  efficient  thermal  shield  materials  are  discussed  in  section  IV-6. 

Here  we  wish  to  present  the  experimental  evaluation  approach  to  the  examina¬ 
tion  of  these  newer  material  formulations.  Included  in  what  is  termed  "newer 
formulations"  can  be  existing  materials  which  look  promising  but  have  not 
received  any  significant  experimental  evaluations. 

As  discussed  earlier, ablation  material  behavior  can  be  characterized 
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aa  steady  state  or  non-steady  state, depending  primarily  upon  the  level  of 
the  heating  rate  and  the  time  duration  of  the  heating  pulse.  Although  both 
areas  are  worthy  of  investigation  in  any  materials  program,  for  general 
screening  of  heat  protection  materials,  evaluations  in  high  intensity,  high 
enthalpy  environments  which  produce  steady  state  response  will  provide  a 
meaningful  yardstick.  This  is  not  to  be  considered  a  downgrading  of  the 
need  for  testing  in  the  non-steady  state  ablation  environments, but  rather  a 
choice  for  the  sake  of  economy  between  the  two  test  environments  in  view 
of  the  extremely  severe  stagnation  point  heating  expected  for  proposed 
missions.  Future  studies  can  further  explore  subsequently  proven  newer 
materials. 

This  study  would  also  be  conducted  in  the  tandem-Gerdien  radiation 
arc  under  similar  conditions  as  imposed  in  the  previously  described  sys¬ 
tems  study.  Criteria  for  the  choice  of  the  test  conditions  would  be  the 
same  as  in  that  work  and  similar  analyses  of  test  data  would  be  performed. 

Choice  of  materials  would  be  based  on  the  outgrowth  of  the  materials 
development  program  discussed  in  Section  IV- 6,  although  this  factor  need 
not  be  very  restrictive.  Sufficient  data  are  currently  available  within  the 
Materials  Section  of  the  SSL  to  establish  early  in  this  program  .during  the 
time  of  test  environment  establishment  and  characterization,  the  types  and 
formulations  of  materials  with  which  to  initiate  testing.  Once  this  initial 
flow  characterization  and  material  testing  have  been  completed  it  is  expect¬ 
ed  that  newly  formulated  materials  emanating  from  the  proposed  materials 
program  would  be  factored  into  this  proposed  experimental  evaluation.  As 
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mentioned  earlier  exieting,  currently  uncharacterized  materials  that  look 
promising  can  also  be  included  in  the  program,  especially  in  the  early 
phases. 

Listed  in  Tables  5. 1  -  5.  5  are  the  nominal  test  conditions  and  the 
number  of  runs  recommended  for  these  studies. 
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TABLE  5. 1  Nominal  Teat  Conditions  Tandem-Gerdien  Radiation  Arc 


‘s/RT 


^RAD, 

BTU/Ft^/Sec. 


BTU/Ft 


CONV 

iFt^/Sec. 


A 

B 

C 

700 

1000 

500 

1000 

3000 

D 

E 

550 

500 

500 

1000 

F 

G 

400 

200 

200 

1000 

TABLE  5.  2  Nominal  Test  Conditions  Tandem -Gerdien  Hypersonic  Arc  Tunnel 

BTU/Ft2,Sec. 


400 


50-75 


TABLE  5.  3  Experimental  Program-Run  Determination  Systems  Study  - 

Radiation  Arc 

1.  Measurement  of  Stagnation  Enthalpy 

3  Test  Points  x  2  ^'^”*/Test  x  3  Gases  18 

2.  Measurement  of  Radiant  Intensities 

3  Test  Points  x  5  ^'^"•/Test  x  3  Gases  45 

3.  Measurement  of  Combined  Radiative  and  Convective  Heat  Transfer 

7  Test  Points  x  3  ^'*”*/TeBt  x  3  Gases  63 

4.  Model  Testing 

Phenolic  -  Nylon 

7  Test  Points  x  6  Testing  Times  x  3  Gases  126 

252 
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TABLE  5.4  Experimental  Program  -  Run  Determination  Syetems  Study 

Hypersonic  Arc  Tunnel 

1.  Measurement  of  Stagnation  Enthalpy 

1  Test  Point  x  2  Runs/ Test  x  3  Gases  6 

2.  Measurement  of  Convective  Heat  Transfer 

1  Test  Point  x  3  Runs/ Test  x  3  Gases  9 

3.  Model  Testing 

Phenolic  -Nylon 

1  Test  Point  x  2  Test  Times  x  3  Runs/ Time 
X  3  Gases  18 

33 


TABLE  5,  5  Experimental  Program  -  Run  Determination  New  Material 

Study  Radiation  Arc 

2  Test  Points  (A,  C) 

X  6  Test  Times  x  3  Gases  x  7  mat'l  252 
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TEST  CONDITIONS 

VrT,-  175  •’vkTo*  175 

Q*  60  BTU/FT*SEC;f»iaOSEC  Q«830  BTU/^T  SEC;  t*IO  SEC 

Qtotal  *  ^200  BTU/fT*  =  8300  BTU/^T* 

Figure  IV-5, 1  Compariion  of  Material  Behavior  at  High  and  Low  Ablation  Rates 
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Figure  IV-5.  2  Stagnation  Point  Heat  Transfer 
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BTU/FT*  SEC 


Figure  IV- 5.  5  Stagnation  Point  Heat  Transfer 
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IV. 6  HEAT  SHIELD  MATERIALS 


6.  1  Introduction 

The  choice  of  thermal  protection  materials  for  Venus-Mars  entry  is 
made  difficult  by  the  extremely  high  heating  rates  and  shear  forces  which 
will  be  encountered.  For  a  90°  entry  into  Venus  with  a  250  lb,  capsule 
(W/CdA  =  60)  the  peak  heating  rates  expected  are  21,700  to  32,  300  BTU/ft.^ 
Preliminary  estimates  (Ref.  18}  indicate  that  as  a  consequence  of  these 
heating  rates,  81  lbs.  of  phenolic-nylon  for  heat  protection  would  be  required 
in  the  250  lb.  vehicle  if  phenolic -nylon  were  the  heat  shield  choice.  With 
the  heat  shield  weight  being  such  a  high  percentage  of  the  total  weight,  any 
reduction  in  weight  will  have  a  strong  influence  on  payload  capability.  Added 
to  the  heating  rate  problem  is  the  fact  that  peak  shear  forces  can  be  expected 
to  reach  60  to  70  lbs.  /ft.  Thus,  the  heat  shield  material  must  have  great 
mechanical  integrity  and  must  not  lose  excess  surface  material  as  a  result 
of  shear  forces  during  the  heat  flux  accompanying  entry,  A  still  further 
complication  is  the  difference  in  Venus  atmosphere  from  that  of  the  Earth, 
i.  e.  ,  greater  percentage  of  CO2  on  Venus. 

These  conditions  are  vastly  different  from  those  for  which  present 
day  heat  shield  materials  have  been  developed.  Therefore  new  approaches 
are  required  which  may  include  not  only  modification  of  present  materials 
but  also  synthesis  of  completely  new  heat  shield  materials  or  composites. 

The  two  requirements  which  appear  to  be  most  prominent  are:  (1)  the  need 
for  a  material  with  an  extremely  high  heat  of  ablation,  and  (2)  the  need  for  a 
reinforced  or  very  strong  material.  The  material  must  also  have  good 
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physical  properties,  and  an  attempt  should  be  made  to  reduce  or  eliminate 
any  physical  incompatibility  between  the  heat  shield  and  the  vehicle  sub¬ 
structure. 

6.  2  Materials  Classes 

Based  on  the  foregoing  requirements  four  classes  of  materials  should 
be  studied  which  either  alone  or  in  combination  appear  to  hold  the  most  pro¬ 
mise  for  success.  These  are: 

(1)  New  organic  resins  with  high  heats  of  ablation. 

(2)  Reinforced  organic  ablation  materials. 

(3)  New  silicone  polymers. 

(4)  Refractory  materials  such  as  pyrolytic  graphite. 

These  classes  of  materials  are  discussed  briefly  in  the  following  paragraphs. 

The  200  Series  material,  which  is  a  quartz  fiber  reinforced  resin, 
has  been  tested  in  a  pebble  bed  heater  jet  (4000°F,  300  BTU/ft^/sec,  900 
psia).  The  weight  loss  was  found  to  be  markedly  less  than  that  of  phenolic 
nylon  and  the  linear  erosion  rate  about  15  times  less  than  that  of  PN.  The 
latter  difference  is  due  partially  to  the  intumescent  character  of  the  Century 
material  as  it  ablates. 

A  second  ablation  condidate  has  shown  a  very  high  heat  of  ablation, 
approximately  three  times  that  of  Series  100  Resin  and  six  times  that  of  PN. 
This  material  involves  the  use  of  a  radically  new  curing  system  with  an  epoxy 
resin.  Such  material  could  be  reinforced  by  high  melting,  high  strength 
fibers.  Modifications  might  include  the  synthesis  of  graft  polymers  using 
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this  material  as  a  backbone.  Prominent  consideration  should  be  given  to 
organometallic  grafts  such  as  silicones. 

A  third  candidate  material  as  a  highly  crosslinked  low  siloxane  resin 
silicone.  Recent  data  indicate  that  this  material  has  a  higher  heat  of  ablation 
than  the  Century  Series  or  PN.  Again  these  materials  could  be  strengthened 
by  reinforcing  agents  and  by  altering  the  molecular  structure. 

A  fourth  possibility  is  the  use  of  pyrolytic -graphite  either  by  itself 
or  in  a  composite.  As  a  first  approximation  it  is  estimated  that  the  weight 
of  PG  required  for  the  vehicle  mentioned  earlier  might  be  one  half  that  for 
PN.  The  insulating  properties  of  PG  coupled  with  the  ablation  properties  of 
materials  mentioned  before  should  result  in  composites  of  great  reliability. 
Such  a  composite  might  take  the  form  of  a  graphite -resin  laminate. 

Vacuum  aging  experiments  with  PN  and  Century  Resins  are  relevant 
to  this  problem  since  a  Venus-Mars  entry  vehicle  will  spend  considerable 
time  in  high  vacuum  before  entry.  These  materials  have  been  exposed  to 
high  vacuum  and  elevated  temperatures  for  extended  periods  of  time.  Under 
these  conditions  no  effect  on  the  heat  of  ablation  of  PN  was  noticed.  However, 
a  slight  improvement  in  the  heat  of  ablation  of  the  Century  material  was  re¬ 
corded. 

The  effect  of  particle  radiation  on  heat  protection  materials  for  a 
Mars -Venus  type  mission  have  not  been  explored.  It  is  possible  that  the 
radiation  encountered  would  have  an  adverse  effect  on  ablation  performance. 

In  going  to  Venus  the  outside  of  the  vehicle  would  receive  10^  to  10^  rads 
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from  solar  emission.  Inside  of  the  first  gm/cm  the  dose  should  be  reduced 
to  10^  to  10^  rads. 

The  above  considerations  should  form  an  integral  part  of  a  materials 
development  program.  The  problems  to  be  encountered  are  such  that  pre¬ 
sent  materials  can  not  be  expected  to  perform  efficiently  and  hence  additional 
research  is  a  necessity. 

6.  3  Materials  Development  and  Evaluation  Program 

A  proposed  program  of  materials  development  would  concentrate  on 
producing  a  heat  shield  material  having  a  high  heat  of  ablation  and  a  high  re¬ 
sistance  to  shear  forces.  Good  physical  properties  would  also  be  emphasized 
and  an  attempt  would  be  made  to  reduce  or  eliminate  any  physical  incompati¬ 
bility  between  the  heat  shield  and  the  vehicle  structure. 

The  program  would  consist  essentially  of  three  tasks  to  run  concurrently 
during  the  course  of  the  study.  These  are:  (1)  Resin  synthesis  and  develop¬ 
ment,  (2)  Formulations  and  reinforcement,  and  (3)  Materials  evaluation  and 
character!  zation. 

In  the  materials  evaluation  and  characterization  portion  of  the  program 
a  complete  range  of  physical  and  thermal  evaluations  will  be  included.  Screen¬ 
ing  tests  will  be  run  to  indicate  the  more  promising  formulations  and  compos¬ 
ites  and  these  will  then  be  tested  further  in  the  arc  facilities  described  in 
previous  sections  of  this  report.  As  this  testing  would  be  conducted  under 
conditions  of  heating  (i.e. ,  radiative/convective)  and  in  atmosphere  pertinent 
to  Venus  entry,  the  results  would  provide  direct  inputs  to  design  calculations 
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of  heat  shield  weights  as  well  as  provide  guidelines  for  further  materials 
improvement. 

The  characterization  of  these  materials  would  be  accomplished  by 
rigorous  examination  of  newly  developed  materials  as  well  as  determination 
of  the  products  of  decomposition  brought  about  by  ablation  testing.  It  is  im¬ 
portant  to  know  the  physical  characteristics  of  the  heat  shield  material  over 
a  wide  temperature  range  in  order  to  determine  design  criteria  as  well  as  to 
determine  physicochemical  interactions  which  might  affect  ablation  perform¬ 
ance.  The  gaseous  and  liquid  decomposition  products  of  the  ablating  material 
determine  the  efficiency  with  which  the  shield  will  perform,  and  must  be 
correlated  with  the  entry  environment  anticipated.  The  structure  of  the  char 
is  also  of  paramount  importance,  since  char  stability  will  determine  the  ex¬ 
tent  to  which  the  plasma  will  erode  the  virgin  plastic  -  char  interface. 

Vacuum  aging  experiments  with  the  more  promising  materials  should 
be  carried  out  since  a  Venus-Mars  entry  vehicle  will  spend  considerable 
time  in  high  vacuum  before  entry.  The  materials  will  be  exposed  to  high 
vacuum  and  elevated  temperatures  for  extended  periods  of  time  to  determine 
the  effect  on  material  degradation  or  on  the  thermal  properties  such  as  that 
of  heat  of  ablation. 

The  foregoing  considerations  of  synthesis,  evaluation  and  characteri¬ 
zation  from  an  integral  part  of  a  materials  development  program.  The  prob¬ 
lems  to  be  encountered  in  Venus-Mars  entry  are  such  that  present  materials 
cannot  be  expected  to  perform  adequately,  and  further  research  is  a  necessity 
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IV. 7  VENUS  IONOSPHERE  PROPAGATION 


A  problem  associated  with  the  determination  of  Venusian  ionospheric 
attenuation  is  that  of  the  effective  electron  collision  cross  section.  Even 
if  an  assumed  model  atmosphere  is  correct  there  is  considerable  uncertainty 
in  the  values  of  the  cross  section.  For  the  expected  electron  energies  well 
below  a  volt,  the  available  data  is  very  sparse  (Ref.  30).  The  uncertainty 
appears  to  be  of  the  order  of  one  magnitude. 

Another  factor  is  that  if  any  new  data  becomes  available  (from  such 
as  a  Venus  Flyby  showing  higher  concentrations  of  CO2  on  Venus),  the 
whole  picture  could  change  and  possible  attenuation  estimates  would  be 
much  higher  and  would  necessitate  an  experimental  program  to  better  deter¬ 
mine  the  effective  collision  cross  section  for  the  Venus  model  atmosphere 
(in  the  energy  range  of  0.001  to  1  ev).  Another  factor  that  should  be  in¬ 
vestigated  is  the  effect  of  ionosphere  electron  concentration  estimate.  The 
present  study  used  an  estimate  of  10^  electrons/cc.  This  is  JPL' s  "best 
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model"  estimate;  their  high  extreme  model  has  an  estimate  of  10  electrons/ 
cc.  This  model  will  make  conditions  much  less  favorable  and  a  short  study 
is  here  indicated  to  determine  what  attenuation  could  be  expected  at  this 
limit. 
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Prediction  of  stagnation  point  heat  transfer  for  Mars 
and  Venus  ballistic  entry  based  on  high  enthalpy 
shock  tube  studies. 

The  development  of  an  experimental  facility,  6  inch 
arc  driven  helium  shock  tube,  capable  of  simulating 
planetary  entry  velocities  up  to  45,000  ft/sec. 

The  performance  of  several  ’blass"  ablation 
materials  relative  to  air  in  simulated  atmospheres  of 
Mars  and  Venus  at  gas  enthalpies  proportional  to 
entry  velocities  of  22,  000  to  28,  000  ft/ sec. 

The  relative  efficiency  of  nitrogen  atoms  and  mole¬ 
cules  in  producing  dissociation  as  determined  by 
studies  of  combustion  generated  shock  structures. 
Development  of  a  versatile  equilibrium  thermo¬ 
chemistry  computer  program  for  planetary  atmos¬ 
pheres. 

Assessment  of  effect  of  Venus  ionosphere  on  commun 
ication  attenuation. 

Recommended  future  experimental  studies. 
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